Partea I Elemente de mecanic¼ a teoretic¼ a 1 Capitolul 1 Cinematica 1 1 1 1 1 Cinematica punctului material Viteza unui punct material Fie O; i1 ; i2 ; i3 un reper cartezian în spa¸tiul tridimensional euclidian E3 : Vom nota cu r = OP = xi1 + yi2 + zi3 (1 1) vectorul de pozi¸tie al punctului generic P: Presupunând c¼ a P î¸si schimb¼ a pozi¸tia în timp, avem r = r(t): (1 2) Rela¸tia (1 2) poate … considerat¼ a ca reprezentând ecua¸tia unei curbe numit¼ a traiectorie În general se presupune c¼ a r(t) are derivate de ordinul doi Viteza punctului P este derivata în raport cu timpul a vectorului de pozi¸tie v(t0 ) = 1 1 2 : dr(t0 ) r(t) = r(t0 ) = lim t!t0 dt t r(t0 ) P (t0 ) P (t) = lim : t!t0 t0 t t0 Accelera¸ tia unui punct material Accelera¸tia punctului material P este derivata vitezei în raport cu timpul, sau echivalent, derivata a doua a vectorului de pozi¸tie în raport cu timpul a(t0 ) = d2 r(t0 ) dv(t0 ) v(t) = r(t0 ) = = lim t!t0 dt2 dt t 3 v(t0 ) : t0 4 CINEMATICA Cum derivata a doua poate … de…nita direct, f¼ ar¼ a calculul prealabil al derivatei întâi, avem de asemenea a(t0 ) = r(t0 ) = lim r (t0 + ) + r (t0 = lim 2r (t0 ) P (t0 ) P (t0 + ) + P (t0 ) P (t0 ) 2 !0 1 1 3 ) 2 !0 = : Proiec¸ tiile vitezei ¸ si accelera¸ tiei pe axele triedrului lui Frenet ata¸ sat traiectoriei Fie s(t) = Z tq x(t)2 + y(t)2 + z(t)2 dt t0 abscisa curbilinie a traiectoriei Viteza punctului material este v(t) = dr(t) dr ds ds = = = s= v dt ds dt dt (1 3) iar accelera¸tia a= dv d2 s d = + dt dt2 ds ds dt 2 v2 = v+ ; (1 4) unde reprezint¼ a tangenta la traiectorie, este normala la traiectorie, este raza de curbur¼ a iar v = kvk este m¼ arimea vitezei 1 1 4 Proiectiile accelera¸ tiei pe axele triedrului lui Darboux În cazul în care punctul material se mi¸sc¼ a pe o suprafa¸ta¼, putem ata¸sa traiectoiei triedrul lui Darboux (!; ! n;! m) ; unde ! este versorul tangentei la traiectorie iar ! n versorul normalei la suprafa¸ta¼ Notând cu unghiul ! ! dintre ¸si n , din (1 4) rezult¼ a 2 2 n m v v ! +! m ; a = !v + ! n unde 1 n = cos curbura geodezic¼ a reprezint¼ a curbura normal¼ a iar (1 5) 1 m = sin repezint¼ a CINEMATICA PUNCTULUI MATERIAL 1 1 5 5 Coordonate curbilinii Presupunem c¼ a putem determina vectorul de pozi¸tie al unui punct în func¸tie de trei parametri independen¸ti (numi¸ti coordonate curbilinii) r = xi1 + yi2 + zi3 = r q 1 ; q 2 ; q 3 astfel încât det @ (x; y; z) = @ (q 1 ; q 2 ; q 3 ) @r @r @r ; ; @q 1 @q 2 @q 3 6= 0: @r @r @r ; ; sunt liniar independen¸ti ¸si putem considera @q 1 @q 2 @q 3 c¼ a determin¼ a o baz¼ a în E3 numit¼ a baza local¼a natural¼a : S¼ a consider¼ am în @r spa¸tiul euclidian E3 reperul local fP; e1 ; e2 ; e3 g cu ei = : Schimbând @q i eventual între ei pe q 1 ¸si q 2 ; putem considera c¼ a acesta este un reper drept A¸sadar vectorii (e1 ; e2 ; e3 ) > 0: Fie în continuare matricea simetric¼ a (gij )1 i;j 3 cu gij = gji = ei ej: Dac¼ a gij = 0 pentru i 6= j spunem c¼ a avem un sistem de coordonate ortogonale Notând cu g = det (gij ) s¼ a demonstr¼ am rela¸tia (e1 ; e2 ; e3 ) = p g: (1 6) Folosind identitatea lui Lagrange avem succesiv (e1 ; e2 ; e3 )2 = [e1 (e2 = e21 (e2 = e21 = h e22 e23 e3 )2 2 (e2 e3 ) i [e1 (e2 [e2 (e1 e3 ) e1 e1 e1 e2 e1 e3 e2 e1 e2 e2 e2 e3 e3 e1 e3 e2 e3 e3 de unde rezult¼ a (1 6) e3 )]2 = e3 )]2 = e3 (e1 e2 )]2 = = det (gij ) = g > 0; 6 1 1 6 CINEMATICA Baze reciproce Componentele contravariante ¸ si covariante ale unui vector din E3 Dac¼ a fe1 ; e2 ; e3 g este o baz¼ a în E3 atunci ¸si e1 ; e2 ; e3 cu e1 = e2 e3 2 e3 e1 3 e1 e2 p ; e = p ; e = p g g g este o baz¼ a, numit¼ a baza reciproc¼a În raport cu baza fe1 ; e2 ; e3 g ; orice vector X 2 E3 se poate scrie sub forma X= 3 X X i ei : i=1 Coordonatele X i poarta numele de componente contravariante al vectorului X: În raport cu baza reciproc¼ a e1 ; e2 ; e3 ; orice vector X 2 E3 se poate scrie sub forma 3 X X= Xi ei : i=1 Coordonatele Xi poart¼ a numele de componente covariante al vectorului X: P În cele ce urmeaz¼a, în aceast¼a sec¸tiune, vom omite s¼a mai folosim semnul ; considerând implicit c¼a dac¼a în cadrul unui produs un indice se repet¼a atunci se face sumare dupa acel indice de la valoarea 1 pân¼a la valoarea 3 a indicelui Fie în continuare matricea simetric¼ a g ij 1 i;j 3 cu g ij = g ji = ei ej : S¼ a ar¼ at¼ am c¼ a g ij = (gij ) 1 ; echivalent cu g ji gik = j k: (1 7) Pentru aceasta este su…cient s¼ a veri…c¼ am c¼ a g ij = Gij g unde Gij reprezint¼ a complementul algebric al lui gij: Vom efectua veri…carea pentru g 12 : Avem succesiv g 12 = e1 e2 = (e2 e3 ) (e3 g e1 ) = (e2 e3 ; e3 ; e1 ) = g CINEMATICA PUNCTULUI MATERIAL 7 (e3 ; e1 ; e2 e3 ) e3 (e1 (e2 e3 )) = = g g g31 g32 g33 g12 G12 e3 [e2 (e1 e3 ) e3 (e1 e2 )] = = = : g g g Vom stabili în continuare regula de coborâre s¸i ridicare a indicelui, adic¼ a regula de calculare a componentelor covariante când cunoa¸stem componentele contravariante ¸si invers Pornim pentru aceasta de la egalitatea u¸sor de veri…cat = ei ej = j i: In rela¸tia X = X j ej = Xk ek (1 8) înmul¸tind scalar cu ei avem X j gij = Xk k i de unde Xi = gij X j : ei Dac¼ a în (1 8) înmul¸tim scalar cu Xj i j (1 9) avem = Xk g ki de unde rezult¼ a X i = g ik Xk : ¸si leg¼ aturile dintre baze pot … exprimate cu ajutorul regulii de ridicare ¸si coborare a indicelui Plec¼ am de la rela¸tia ei = j je : Înmul¸tind scalar cu ek avem gik = j j k = k deci ei = gik ek : Înmul¸tind în (1 10) cu g ji (1 10) (cu sumare subîn¸teleas¼ a dup¼ a i) avem g ji ei = g ji gik ek de unde, în virtutea lui (1 7) rezult¼ a ej = g ji ei : (1 11) Observa¸tie Dac¼ a fP; e1 ; e2 ; e3 g este un reper cartezian, atunci gij = i = j iar componentele covariante ¸si contravariante ale unui vector coincid ¸si acela¸si lucru se întâmpl¼ a cu baza ini¸tial¼ a ¸si baza reciproc¼ a g ij 8 CINEMATICA 1 1 7 Componentele …zice ale unui vector Prin normarea bazei X = X j ej = 3 X j=1 X p ej ek p = Xk ek = Xk g kk p X j gjj p gjj g kk k=1 se asociaz¼ a …ec¼ arei componente a unui vector o m¼ arime numit¼ a componenta …zic¼ a: p p X j(f iz) = X j gjj ; Xj(f iz) = Xk g kk (f ara sumare) 1 1 8 Componentele contravariante ¸ si covariante ale vitezei ¸ si accelera¸ tiei Din rela¸tia v = v i ei = i @r i q = q ei i @q calcul¼ am componentele contravariante ale vitezei în cazul sistemului de coordonate curbilinii q 1 ; q 2 ; q 3 : i v i = q ; i = 1; 2; 3: (1 12) Cu regula de coborâre a indicelui rezult¼ a componentele covariante ale vitezei j vi = gij q ; i = 1; 2; 3: (1 13) Pentru a calcula componentele accelera¸tiei ¸tinem cont de rela¸tia a= d dv = dt dt @r j q @q j = @2r j k @r j q q + q : @q j @q j @q k Pentru a calcula componentele covariante ale accelera¸tiei, în rela¸tia a = aj ej ; înmultind scalar cu ei ob¸tinem j k j ai = a ei = q ej ei + q q ei echivalent cu j j k ai = q gji + q q i;jk ; @ej @q k i = 1; 2; 3; (1 14) CINEMATICA PUNCTULUI MATERIAL 9 unde i;jk @ej @q k = ei reprezint¼ a coordonatele covariante ale lui (1 15) @ej ¸si poart¼ a numele de simbolul @q k lui Christo¤ el de spe¸ta întâi Observ¼ am c¼ a i;jk = ei @ej @2r @ek = e = ei = i k j k @q j @q @q @q i;kj : (1 16) Pentru a calcula componentele contravariante ale accelera¸tiei folosim regula de ridicare a indicelui: al = g li ai = q j l j j k + q q g li i;jk : (1 17) Din egalitatea @gij @ej @ei = k ej + k ei k @q @q @q ¸si din (1 15) deducem pentru simbolul lui Christo¤el de spe¸ta întâi formula de calcul @gjk @gij 1 @gik + : + k;ij = j k 2 @q @q i @q Introducând simboliul lui Christo¤ el de spe¸ta a doua: l jk = g li i;jk = g li ei @ej @ej = el k @q @q k avem alt¼ a expresie pentru componentele contravariante ale accelera¸tiei: l j k l jk : al = q + q q (1 18) Vom expune în continuare o noua metod¼ a de calcul a componentelor covariante ale accelera¸tiei Introducând energia punctului de mas¼ a egal¼ a cu unitatea j i i 1 1 j 1 T = v2 = q ej q ei = gji q q ; 2 2 2 s¼ a demonstr¼ am c¼ a " # d @T aj = dt @ q: j @T : @q j (1 19) 10 CINEMATICA Avem succesiv # " d @T dt @ q: j @T d :i = gji q j @q dt = @gji : k : i :: i q q + gji q k @q :: i = gji q + i;jk + j;ik j;ik + :: i = gji q + 1 @ :k :i (ei ek ) q q = j 2 @q :: i = gji q + 1 2 1 2 = 1 @gki : k : i q q = 2 @q j @ (ei ej ) @q k :: i = gji q + k i 1 gki q q 2 @ @q j 1 2 i;kj i;jk :k :i j;ik q q 1 2 :k :i k;ij q q = :k :i k;ji q q = = aj : Aplica¸tii 1) S¼a se determine componentele vitezei s¸i accelera¸tiei în coordonate cilindrice Rezolvare Fie coordonatele cilindrice ( ; ; z) ; 2 (0; 1) ; 2 [0; 2 ) ; z 2 R introduse prin rela¸tiile x = cos ; y = sin ; z = z: Fie vectorul de pozi¸tie r = xi1 +yi2 +zi3 ¸si baza local¼ a natural¼ a fe1 ; e2 ; e3 g = fe ; e ; ez g introdus¼ a prin rela¸tiile e = @r @r = cos i1 + sin i2 ; e = = @ @ sin i1 + cos i2 ; ez = i3 : Coordonatele cilindrice sunt ortogonale deoarece 0 1 0 1 1 0 0 (gij )1 i;j 3 = @ 0 2 0 A ; : g ij 1 i;j 3 = @ 0 0 0 1 0 0 1 2 0 1 0 0 A: 1 Componentele contravariante ale vitezei sunt date de rela¸tia : : : v = e + e + zez iar componentele contravariante …zice ale vitezei le sunt date de rela¸tia : v= i + : : i + ziz CINEMATICA PUNCTULUI MATERIAL cu i =e ; i = e 11 ; iz = ez : Componentele covariante ale vitezei sunt date de rela¸tia : 2 v= e + : : e + zez iar componentele covariante …zice ale vitezei le sunt date de rela¸tia : : v= e + : : : + zez = i + e : i + ziz ¸si coincid cu componentele contravariante …zice Fie energia cinetic¼ a a punctului material cu masa egal¼ a cu unitatea 1 1 T = v2 = 2 2 :2 2 + :2 :2 +z : Componentele covariante ale accelera¸tiei sunt a = a = d dt d dt @T : @ @T = @ @T @T = @ : @ d az = dt @T : @z :2 :: 2 :: +2 ; : : ; @T :: = z; @z iar cele contravariante a =a ; a = 2 a ; az = az : Componentele …zice ale accelera¸tiei se ob¸tin din rela¸tia a= :: :2 :: i + +2 :: :: i + ziz : Observa¸tie În cazul în care mi¸scarea punctului material are loc într-un plan (s¼ a zicem planul Oxy ) avem = r = krk ; z = 0 iar componentele …zice ale vitezei ¸si accelera¸tiei se ob¸tin din rela¸tiile : : v = ri + r i ; 12 CINEMATICA :: a= r r :2 :: :: i + r + 2r i : Coordonatele (r; ) poart¼ a numele de coordonate polare 1) S¼a se determine componentele vitezei s¸i accelera¸tiei în coordonate sferice Rezolvare Fie (r; ; ') ; r 2 (0; 1) ; 2 ; ' 2 [0; 2 ), coordonatele sferice introduse prin rela¸tiile x = r sin cos '; y = r sin sin '; z = r cos : Baza local¼ a natural¼ a este format¼ a din vectorii e1 = er = sin cos 'i1 + sin sin 'i2 + cos i3 ; e2 = e = r cos cos 'i1 + cos sin 'i2 e3 = e' = r sin sin i3 ; sin 'i1 + cos 'i2 : Avem de asemenea 1 0 1 0 0 A ; g ij 0 (gij )1 i;j 3 = @ 0 r2 2 2 0 0 r sin 1 i;j 3 0 1 @ 0 = 0 0 1 r2 0 0 0 1 r2 sin 2 1 A; coordonatele sferice …ind deci coordonate ortogonale Vectorii bazei reciproce sunt er = sin cos 'i1 + sin sin 'i2 + cos i3 ; 1 cos cos 'i1 + cos sin 'i2 sin i3 ; r 1 e' = sin 'i1 + cos 'i2 : r sin Componentele contravariante ¸si componentele contravariante …zice ale vitezei se ob¸tin din rela¸tiile : e = : : : : : v = rer + e + 'e' = rer + r e r : + r sin ' e' r sin : Componentele covariante ¸si componentele covariante …zice ale vitezei se ob¸tin din rela¸tiile : : : : : v = rer + r2 e + r2 sin 2 'e' = rer + r : re : + r sin ' (r sin e' ) : CINEMATICA SOLIDULUI RIGID 13 Energia cinetic¼ a a punctului cu masa egal¼ a cu unitatea este 1 1 T = v2 = 2 2 :2 r + r2 :2 :2 + r2 sin 2 ' iar coordonatele covariante ale accelera¸tiei sunt ar = a = d dt d dt 1 2 r :2 :: @T :: = r2 + 2rr @ @T : @ a' = @T :: =r @r @T : @r d dt @T : @' :2 r sin 2 ' ; :2 r2 ' sin cos ; d 2 @T : = r sin 2 ' : @' dt Cinematica solidului rigid 1 2 1 Vectorul de rota¸ tie instantanee Formulele lui Poisson Fie în spa¸tiul euclidian E3 reperul cartezian …x O1 ; i1 ; j 1 ; k 1 : Fie de asemenea reperul mobil O; i; j; k solidar legat de corpul solid rigid (S) (vezi …gura 1 1) Notând cu M un punct al solidului, s¼ a consider¼ am vectorii de pozi¸tie ! ! ! ! r1 = O1 M ; ! r = OM ; r! O = O1 O; astfel ca r1 = x1 i1 + y1 j 1 + z1 k 1 ; r = xi + yj + zk; rO = xO i + yO j + zO k: În timpul mi¸sc¼ arii solidului coordonatele x; y; z sunt constante iar coordonatele xO ; yO ; zO precum ¸si versorii i; j; k variaz¼ a cu timpul Consider¼ am rela¸tiile di = ! 11 i + ! 12 j + ! 13 k; dt (1 20) dj = ! 21 i + ! 22 j + ! 23 k; dt (1 21) dk = ! 31 i + ! 32 j + ! 33 k: dt (1 22) 14 CINEMATICA Figure 1 1: Mi¸scarea solidului rigid ¸tinând cont de condi¸tiile de ortonormalitate 2 2 2 i = 1; j = 1; k = 1; i j = 0; j k = 0; k i = 0 deducem i di dj dk = 0; j = 0; k = 0; dt dt dt (1 23) dj di dk dj di dk +j = 0; j +k = 0; k +i = 0: dt dt dt dt dt dt Din (1 20) - (1 24) rezult¼ a (1 24) i ! 11 = ! 22 = ! 33 = 0; ! 12 + ! 21 = 0; ! 13 + ! 31 = 0; ! 23 + ! 32 = 0: Am de…nit în felul acesta tensorul antisimetric ! ! = (! ij )1 ¸si tensorul rota¸tie Punând dj k = ! 23 = p(t); dt i dk = ! 31 = q(t); dt j di = ! 12 = r(t) dt i;j 3 numit CINEMATICA SOLIDULUI RIGID 15 se de…ne¸ste un vector ! = pi + qj + rk numit vectorul rota¸tie (vitez¼a unghiular¼a) Cu noile nota¸tii avem di = rj dt qk = ! i; (1 25) dj = pk dt ri = ! j; (1 26) dk = qi dt pj = ! k: (1 27) Formulele (1 25) - (1 27) poart¼ a numele de formulele lui Poisson 1 2 2 Distribu¸ tia vitezei ¸ si accelera¸ tiei în cazul mi¸ sc¼ arii solidului rigid Vom demonstra mai întâi rela¸tia dr =! dt r: (1 28) Avem succesiv dr di dj dk =x +y +z = x! dt dt dt dt i + y! j + z! k=! r: Viteza punctului generic M al solidului, raportat¼ a la reperul …x va … v(t) = dr1 drO dr = + ; dt dt dt de unde, notând cu v O viteza punctului O avem v = vO + ! r: (1 29) Accelera¸tia punctului M este a= dv dv O d! = + dt dt dt r+! dr ; dt sau notând cu aO accelera¸tia punctului O rezult¼ a : a = aO + ! r+! (! r) : (1 30) 16 CINEMATICA 1 3 1 3 1 Teoria mi¸ sc¼ arii relative Cazul unui reper …x ¸ si al unui reper mobil Vom considera reperul …x O1 ; i1 ; j 1 ; k 1 ¸si reperul mobil O; i; j; k : Notând cu M un punct oarecare vom folosi nota¸tiile r1 = O1 M = x1 i1 + y1 j 1 + z1 k 1 ; r = OM = xi + yj + zk; rO = O1 O: Vom de…ni viteza relativ¼ a (adic¼ a viteza în raport cu reperul mobil O; i; j; k ) a punctului M prin rela¸tia vr = @r : : : = xi + yj + zk; @t (1 31) @ indic¼ am o derivare par¸tial¼ a în sensul ca vectorii i; j; k nu sunt @t deriva¸ti ) Viteza absolut¼ a (adica viteza în raport cu reperul …x O1 ; i1 ; j 1 ; k 1 ) a punctului M va … (prin va = drO dr dr1 d = + = vO + xi + yj + zk = dt dt dt dt : : : : : : = v O + xi + y j + z k + xi + yj + zk sau, ¸tinând cont de (1 31) ¸si de formulele lui Poisson, va = vO + ! r + vr (1 32) unde ! este vectorul ro¸tatie al reperului mobil O; i; j; k în raport cu reperul …x O1 ; i1 ; j 1 ; k 1 : Introducând viteza de transport (ca …ind viteza unui punct solidar legat de reperul mobil): vt = vO + ! r putem rescrie rela¸tia (1 32) sub forma va = vt + vr : (1 33) Pentru a ob¸tine expresia accelera¸tiei absolute aa a punctului M , deriv¼ am în raport cu timpul rela¸tia (1 32): aa = : dv a dv O = +! dt dt r+! dr dv r + : dt dt (1 34) ¼ TEORIA MI¸ SCARII RELATIVE 17 Vom introduce accelera¸tia relativ¼ a (adic¼ a accelera¸tia în raport cu reperul mobil O; i; j; k ) a punctului M prin rela¸tia ar = @v r :: :: :: = xi + yj + zk: @t (1 35) Tinând ¸ cont de (1 34) ¸si de formulele lui Poisson ob¸tinem dv r d : : : xi + yj + zk = ! = dt dt Cum dr =! dt v r + ar : (1 36) r + vr ; din (1 34) - (1 36) rezult¼ a : aa = aO + ! r+! (! unde aO = r) + 2! v r + ar ; (1 37) dv O dt reprezint¼ a accelera¸tia punctului O: Introducând accelera¸tia de transport (ca …ind accelera¸tia unui punct solidar legat de reperul mobil): : at = aO + ! r+! (! r) ; ¸si accelera¸tia complementar a¼ (sau accelera¸tia lui Coriolis) ac = 2! vr ; putem rescrie relatia (1 37) sub forma aa = at + ac + ar : (1 38) 18 CINEMATICA Capitolul 2 Dinamica punctului ¸ si a sistemelor de puncte materiale 2 1 Legile mecanicii clasice Mecanica newtonian¼ a clasic¼ a admite ca realitate …zic¼ a existen¸ta unui spa¸tiu absolut ¸si a unui timp absolut , independente între ele ¸si f¼ ar¼ a nici-o rela¸tie cu materia, în sensul c¼ a modelul matematic al acestui spa¸tiu este spa¸tiul euclidian R3 , iar oric¼ arui eveniment mecanic, care se traduce prin prezen¸ta unei particule materiale într-un punct al acestui spa¸tiu la un moment dat, îi corespunde un element al produsului cartezian R3 R: În mecanica clasic¼ a se postuleaz¼ a existen¸ta unui sistem de referin¸ta¼ numit reper iner¸tial astfel ca, în cadrul geometriei euclidiene s¼ a …e veri…cate urm¼ atoarele principii enun¸tate de Newton: 1 Principiul iner¸ tiei Acest principiu a…rm¼ a c¼ a dac¼a asupra unui punct material nu actioneaz¼a nici-o for¸t¼a, el se mi¸ sca rectiliniu s¸i uniform sau este în stare de repaus fa¸t¼a de reperul iner¸tial 2 Legea a doua a lui Newton Conform acestei legi for¸ta F care ac¸tioneaz¼a asupra unui punct material este egal¼a cu produsul dintre masa m a punctului s¸i accelera¸tia sa (notat¼a cu a) : F = ma (2 1) 3 Principiul ac¸ tiunii ¸ si reac¸ tiunii Acest principiu postuleaz¼ a c¼ a dac¼a ac¸tiunea exercitat¼a de punctul material M2 asupra lui M1 se 19 20 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE reprezint¼a prin for¸ta F12 s¸i dac¼a F21 este for¸ta exercitat¼a de punctul material M1 asupra lui M2 , atunci F12 s¸i F21 au ca suport dreapta M1 M2 s¸i F12 + F21 = 0: (2 2) 2 1 1 Principiul de relativitate în mecanica clasic¼ a Ecua¸tia (2 1) este valabil¼ a fa¸ta¼ de un reper iner¸tial O1 x1 y1 z1 a c¼ arui existen¸ta¼ a fost postulat¼ a S¼ a consider¼ am un nou reper Oxyz animat fa¸ta¼ de O1 x1 y1 z1 de o mi¸scare de transla¸tie rectilinie ¸si uniform¼ a de vitez¼ a (evident constant¼ a) vO : Cum vectorul de rota¸tie al reperului Oxyz în raport cu reperul O1 x1 y1 z1 este nul, conform legilor de compunere a vitezelor ¸si accelera¸tiilor, între viteza v ¸si accelera¸tia a a unui punct M fa¸ta¼ de reperul O1 x1 y1 z1 ¸si viteza relativ¼ a vr ¸si accelera¸tia relativ¼ a ar fa¸ta¼ de noul reper, avem rela¸tiile v = vr + vO (2 3) a = ar (2 4) Din rela¸tia (2 3) deducem c¼ a principiul iner¸tiei este valabil ¸si fa¸ta de reperul Oxyz iar din rela¸tia (2 4) rezult¼ a c¼ a fa¸ta¼ de acest reper isi p¼ astreaz¼ a valabilitatea ¸si legea a doua a lui Newton c¼ aci F = mar : (2 5) Rezulta de aici c¼ a dac¼a avem un reper iner¸tial, exist¼a o in…nitate de astfel de repere iner¸tiale, deducându-se unul din altul printr-o transla¸tie rectilinie s¸i uniform¼a Aceast¼ a a…rma¸tie constituie principiul de relativitate al lui Galileu 2 2 Mi¸ scarea punctului material liber Legea a doua a lui Newton, scris¼ a pe coordonate sub forma m• x = F1 (x; y; z; x; y; z; t); m• y = F2 (x; y; z; x; y; z; t); m• z = F3 (x; y; z; x; y; z; t) constituie un sistem de 3 ecua¸tii diferen¸tiale ordinare de ordinul doi (2 6) MI¸ SCAREA PUNCTULUI MATERIAL LIBER 2 2 1 21 Legi de conservare în dinamica punctului material Integrarea sistemului (2 6) este u¸surat¼ a de cunoasterea unei sau mai multor integrale prime, adica rela¸tii de forma f (x; y; z; x; y; z; t) = C (2 7) în care C = const: iar (x; y; z) reprezint¼ a o solu¸tie a sistemului (2 6) Putem avea cel mult 6 integrale prime independente de forma fi (x; y; z; x; y; z; t) = Ci ; i = 1; 2; :::; 6 (2 8) Integralele prime sunt independente dac¼ a avem rela¸tia det @(f1 ; f2 ; f3 ; f4 ; f5 ; f6 ) 6= 0: @(x; y; z; x; y; z) (2 9) În acest caz, cu teorema func¸tiilor implicite, din (2 9) se ob¸tine x = x(t; C1 ; :::; C6 ); y = y(t; C1 ; :::; C6 ); z = z(t; C1 ; :::; C6 ) (2 10) x = x(t; C1 ; :::; C6 ); y = y(t; C1 ; :::; C6 ); z = z(t; C1 ; :::; C6 ) (2 11) constantele C1 ; :::; C6 ob¸tinându-se din condi¸tiile ini¸tiale Integralele prime ale sistemului (2 6), con¸tinând rela¸tii între coordonatele punctului material, componentele vitezei ¸si timp, poart¼ a numele de legi de conservare 2 2 2 Teoremele generale ale mecanicii punctului material Teoremele generale ale mecanicii punctului material sunt teorema impulsului, teorema momentului cinetic ¸si teorema energiei cinetice Ele sunt importante printre altele pentru c¼ a în anumite cazuri ne pot furniza integrale prime ale sistemului (2 6) Teorema impulsului Deoarece masa punctului material este constant¼ a se deduce ma = d(mv) = F: dt Introducând impulsul (sau cantitatea de mi¸ scare) ca produsul dintre masa punctului material ¸si viteza sa H = mv ob¸tinem: (2 12) 22 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE 1 Teorema impulsului Derivata impulsului în raport cu timpul este egal¼a cu for¸ta : dH =F (2 13) dt Dac¼ a for¸ta F îndepline¸ste anumite condi¸tii, putem deduce de aici integrale prime ale sistemului (2 6) De exemplu, daca for¸ta F este perpendicular¼ a pe un vector …x u( ; ; ) vom avea din (2 13): d (H u) = F u = 0 ) H u = C; dt C …ind o constant¼ a Proiec¸tiile vectorului H pe axe …ind Hx = mx; Hy = my; Hz = mz integrala prim¼ a ob¸tinut¼ a se scrie m( x + y + z) = C: Teorema momentului cinetic Fie r =OM vectorul de pozitie al punctului material M în raport cu reperul iner¸tial Oxyz Avem a¸sadar v= dr d2 r ; a= 2: dt dt Deducem de aici mr d d2 r = m(r 2 dt dt dr )=r dt F: Vectorul K=r mv = r H este momentul lui H, aplicat în M fa¸ta de originea O El se nume¸ste moment cinetic sau moment al cantit¼at¸ii de mi¸ scare fa¸t¼a de punctul O Pe de alt¼ a parte, r F este momentul MO al for¸tei F , aplicate în punctul M , fa¸ta¼ de punctul O Avem deci Teorema momentului cinetic Derivata momentului cinetic fa¸t¼a de un punct …x O este egal¼a cu momentul for¸tei fa¸t¼a de acel punct, MI¸ SCAREA PUNCTULUI MATERIAL LIBER 23 dK (2 14) = MO : dt Din considera¸tiile f¼ acute mai sus, observ¼ am c¼ a teorema momentului cinetic este valabil¼ a nu numai fa¸ta¼ de originea O a reperului, ci ¸si fa¸ta¼ de orice punct …x în raport cu reperul iner¸tial considerat Dac¼ a îns¼ a se iau momentele fa¸ta¼ de un punct mobil, enun¸tul teoremei se modi…c¼ a Fie punctul mobil A de vector de pozi¸tie rA (t) ¸si viteza vA : Notând prin KO momentul lui mv fa¸ta¼ de O ¸si prin KA momentul lui mv fa¸ta¼ de A, apoi cu MO momentul lui F fa¸ta¼ de O ¸si cu MA momentul s¼ au fa¸ta¼ de A; avem KO = KA + OA MO = MA + OA mv = KA + rA H; F = MA + rA F: Prin urmare rela¸tia (2 14) conduce la dKA + vA dt H + rA dH = MA + rA dt F: ¸tinând cont de (2 13) deducem dKA = MA dt vA H aceast¼a rela¸tie reprezentând teorema momentului cinetic fa¸t¼a de un punct mobil A: Teorema ariilor Daca momentul MO este perpendicular pe versorul …x u( ; ; ) atunci din (2 14) se deduce d (KO u) = MO u = 0 dt deci KO u = mC unde C este o constanta Avem o integral¼ a prim¼ a Fie Kx = m(y z z y); Ky = m(z x xz); Kz = m(xy y x) proiec¸tiile vectorului KO pe axele de coordonate Integrala prim¼ a ob¸tinut¼ a se mai scrie Kx + Ky + Kz = mC: 24 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE Considerând c¼ a u este versorul axei Oz, vom avea integrala prim¼ a se va scrie mai simplu (xy = y x) = C: = 0; = 1 ¸si (2 15) Fie M (x; y; z) punctul mobil considerat, de vitez¼ a v ¸si P (x; y; 0), proiec¸tia sa pe planul Oxy, care este perpendicular pe versorul u: Punctul P va descrie în planul Oxy o traiectorie care este proiec¸tia ortogonal¼ a pe acest plan a traiectoriei punctului M ; viteza lui P este vP (x; y; 0), proiec¸tia ortogonal¼ a a lui v: Fie Sz (t) suprafa¸ta descris¼ a de raza vectoare OP de la momentul ini¸tial t0 , când M era în M0 iar P era în P0 pân¼ a la momentul t: Not¼ am cu Az (t) aria suprafe¸tei Sz (t) ¸si cu @Sz (t) frontiera ei Pornind de la rela¸tia Z Z Az (t) = dxdy Sz (t) ¸si folosind formula lui Stokes Z Z @Q @R @x @y Sz (t) dxdy = x ;R= 2 Q= 1 Az (t) = 2 Rdx + Qdy @Sz (t) cu ob¸tinem Z Z y 2 xdy ydx: @Sz (t) Cum @Sz (t) este alc¼ atuit din segmentele OP; P0 O ¸si din arcul P0 P , ob¸tinem în de…nitiv, Z 1 Az (t) = xdy ydx (2 16) 2 P0 P deoarece avem dy y = ; (x; y) 2 OP [ P0 O: dx x Din (2 16) rezult¼ a 1 Az (t) = 2 Z t t0 (xy y x)dt: MI¸ SCAREA PUNCTULUI MATERIAL LIBER 25 ¸tinând cont de (2 15) ob¸tinem Az = C (t 2 t0 ): (2 17) Prin urmare avem Teorema ariilor Dac¼a momentul for¸tei F este mereu nul fa¸t¼a de axa Oz, atunci în proiec¸tie pe planul Oxy, aria descris¼a de raza vectoare OP , de la pozi¸tia sa ini¸tial¼a OP0 ; este proportional¼a cu intervalul de timp parcurs dAz 1 M¼ arimea = (xy y x) poart¼ a numele de viteza areolar¼a a punctului dt 2 P Dac¼ a momentul for¸tei este mereu nul fa¸ta¼ de o alt¼ a dreapt¼ a trecând prin origine, de direc¸tie …x¼ a, teorema ariilor se aplic¼ a în proiec¸tie pe un plan perpendicular pe dreapt¼ a ¸si care trece prin origine Constanta C se numeste constanta ariilor Teorema energiei cinetice S¼ a relu¼ am ecua¸tia de mi¸scare (2 1) sub forma dv =F (2 18) dt Înmul¸tind ambii membrii scalar cu v ¸si integrând între momentele t0 (când viteza ia valoarea v0 ) ¸si t, ob¸tinem Z t 1 1 2 2 mv mv = F vdt (2 19) 2 2 0 t0 m M¼ arimea 1 1 T = mv 2 = mv 2 2 2 poart¼ a numele de energie cinetica arimea F v R t a punctului material , m¼ de putere a for¸tei iar m¼ arimea t0 F vdt de lucru mecanic total al for¸tei Conform celor de mai sus avem Teorema energiei cinetice Variatia energiei cinetice între momentele t0 s¸i t este egal¼a cu lucrul mecanic total al for¸tei între aceste momente Putem formula teorema energiei cinetice ¸si sub forma diferen¸tial¼ a dac¼ a în rela¸tia (2 18) înmul¸tim cu dr = vdt: Ob¸tinem atunci rela¸tia d 1 mv 2 2 = F dr (2 20) care ne arat¼ a c¼ a diferen¸tiala energiei cinetice a punctului material este egal¼ a cu lucrul mecanic elementar al for¸tei F dr: 26 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE Câmp de for¸ te care admite un poten¸ tial Dac¼ a for¸ta F nu depinde decât de r, adic¼ a de pozi¸tia punctului material, avem F = F (x; y; z): Se spune c¼ a în acest caz avem de-a face cu un câmp de for¸te În cazul în care câmpul de for¸te admite un poten¸tial U (x; y; z), adic¼ a F = gradU (x; y; z) rela¸tia (2 19) devine 1 mv 2 2 1 mv 2 = 2 0 Z t t0 ( @U dx @U dy @U dz + + )dt = @x dt @y dt @z dt Z t t0 dU dt dt echivalent¼ a cu 1 mv 2 2 1 U (x; y; z) = mv02 2 U (x0 ; y0 ; z0 ) (2 21) Se obi¸snuie¸ste ca func¸tia V = U s¼ a …e denumit¼ a energia poten¸tial¼a a punctului material iar func¸tia T +V energia mecanic¼ a total¼ a A¸sadar rela¸tia (2 21) se mai scrie T + V = const: (2 22) Avem deci legea de conservare : În cazul unui câmp de for¸te care admite un poten¸tial, energia mecanic¼a total¼a a punctului material supus ac¸tiunii acelui câmp este constant¼a Câmpul de for¸te se nume¸ste din acest motiv câmp de for¸te conservativ 2 2 3 Mi¸ scarea unui punct material sub ac¸ tiunea unei for¸ te centrale Spunem c¼ a o for¸ta este central¼a dac¼ a suportul s¼ au trece printr-un punct …x numit pol sau centru Propozi¸ tie Traiectoria unui punct material care se mi¸ sc¼a sub ac¸tiunea unei for¸te centrale este plan¼a Demonstra¸tie Fie O polul for¸tei Din rela¸tia : K O = MO = 0 deducem KO = r mv = C = const: MI¸ SCAREA PUNCTULUI MATERIAL LIBER 27 A¸sadar C r = 0 deci mi¸scarea are loc într-un plan perpendicular pe vectorul C: Considerând reperul Oxyz = O; i1 ; i2 ; i3 astfel încât C k i3 , rezult¼ a c¼ a mi¸scarea are loc în planul Oxy: Trecând la coordonate polare, x = r cos ; y = r sin deducem din teorema ariilor c¼ a : 2 2 4 : : y x = r2 = C: xy Formula lui Binet În cazul unei for¸te centrale avem a¸sadar F = Fi : (2 23) ¸tinând cont de expresiile coordonatelor …zice ale accelera¸tie în coordonate polare din (2 23) deducem în virtutea legii a doua a lui Newton :: m r :: :: m 2r + r =0) r :2 = F; : d r2 dt (2 24) : = 0 ) r2 = C: (2 25) : : Observ¼ am c¼ a reg¼ asim teorema ariilor cu C = r(t0 )2 (t0 ) = r02 0 : Din (2 24) ¸si (2 25) deducem rela¸tia :: mr = F + C2 : r3 (2 26) Cautând ecua¸tia traiectoriei plane sub forma r = r( ); avem : dr : C dr r= = 2 = d r d : :: r= C d2 d 2 1 r : = C d d 1 r C 2 d2 r2 d 2 ; 1 r (2 27) : (2 28) Din (2 26) - (2 28) deducem formula lui Binet: F = mC 2 1 d2 + r2 r d 2 1 r : (2 29) 28 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE Dac¼ a F = F (r; ) cu formula lui Binet se poate a‡a traiectoria punctului material care se mi¸sc¼ a sub actiunea for¸tei centrale Urmeaza apoi s¼ a se a‡e = (t); adica legea orar¼a a mi¸sc¼ arii pe traiectorie Aplica¸tie Un punct material M de mas¼a m este atras de centrul O cu km o for¸t¼a F = unde r = OM iar k este o constant¼a cunoscut¼a Ce r5 vitez¼a ini¸tial¼a v 0 paralel¼a cu axa Oy trebuie s¼a se imprime punctului M a‡at în A 2 Ox pentru ca el s¼a descrie un cerc de diametru OA (se va exprima în func¸tie de k s¸i R)? Dup¼a cât timp ajunge mobilul în O ? Rezolvare Deoarece punctul M descrie semicercul de raz¼ a R (…g 2 1) avem : v0 r = 2R cos ; C = r2 = (2R)2 = 2Rv0; 2R de unde 1 d 1 = ; r 2R cos d 1 r = sin d2 ; 2 2R cos d 2 1 r = 1 + sin2 : 2R cos3 Înlocuind în formula lui Binet ob¸tinem 1 + sin2 2R cos3 + 1 2R cos = k k 2 3 2 ) v0 = 32R2 : (2R cos ) (2Rv0 ) Figure 2 1: Aplica¸tie la formula lui Binet Timpul dup¼ a care mobilul ajunge în O se calculeaz¼ a din legea ariilor : C = r2 = 2 A A R2 R )t=2 = = : t C 2Rv0 2v0 ¼ PUNCTUL MATERIAL SUPUS LA LEGATURI 29 2 3 Mi¸ scarea punctului material supus la leg¼ aturi 2 3 1 Cazul leg¼ aturilor ideale (f¼ ar¼ a frecare) Mi¸ scarea pe o suprafa¸ ta ¼ În cazul mi¸sc¼ arii unui punct material pe suprafa¸ta ' (x; y; z; t) = 0, din legea a doua a lui Newton ob¸tinem ma = F + N (2 30) unde F este rezultanta for¸telor active iar for¸ta de leg¼ atur¼ a (reac¸tiune) este N = grad': În cazul unei suprafe¸te …xe de ecua¸tie '(x; y; z) = 0; notând cu r = r(s); ecua¸tia traiectoriei (unde s reprezint¼ a abscisa curbilinie), se consider¼ a dr grad' triedrul lui Darboux ( ; m; n) cu = ; n = ; m = n : ds kgrad'k Cunoscând proiec¸tiile accelera¸tiei pe axele triedrului lui Darboux a= v+n v2 +m n v2 m ¸si considerând rezultanta for¸telor active F =F + Fm m + Fn n; ob¸tinem ecua¸tiile de mi¸scare în proiec¸tie pe axele triedrului lui Darboux : mv = F ; m v2 = Fn + N; (2 31) (2 32) n m v2 = Fm: (2 33) m S¼ a presupunem în continuare c¼ a for¸ta F este conservativ¼ a având poten¸tialul U (x1 ; x2 ; x3 ) (F = grad U ) Înmul¸tind în (2 31) cu v avem dT @U dx1 @U dx2 @U dx3 dU : = mv v = vF = F v = + + = ; (2 34) dt @x1 dt @x2 dt @x3 dt dt de unde rezult¼ a integrala prim¼ a a energiei T U = const: (2 35) 30 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE Mi¸ scarea pe o suprafa¸ ta a dat¼ a sub form¼ a parametric¼ a ¼ …x¼ Fie r = r q 1 ; q 2 ecua¸tia parametric¼ a a suprafe¸tei Avem rela¸tiile : :1 :2 r = r0q1 q + r0q2 q ; :: :: 1 :: 2 :1 :1 :2 :2 r = r0q1 q + r0q2 q + r00q1 q1 (q )2 + 2r00q1 q2 q q + rq002 q2 (q )2 : (2 36) Tinând ¸ cont de rela¸tiile (2 36), avem în virtutea rela¸tiilor (??) din anex¼ a (formulele lui Gauss) :: :: 1 :: 2 r r0qi = g1i q + g2i q + :1 2 i;11 (q ) :: +2 :1 :2 i;12 q q :1 :1 :2 + :2 2 i;22 (q ) ; i = 1; 2: (2 37) :2 r U = G11 (q )2 + 2G12 q q + G22 (q )2 : (2 38) Cu acestea ecua¸tia (2 30) în proiec¸tie pe suporturile tangentelor r0q1 ; r0q2 r0q1 r0q2 la suprafa¸ta¼ ne d¼ a ¸si normalei U = r0q1 r0q2 :: 1 :: 2 m g1i q + g2i q + :1 2 i;11 (q ) +2 :1 :2 i;12 q q + :2 2 i;22 (q ) = Qi ; i = 1; 2 (2 39) cu Qi = F r0qi ¸si :1 :1 :2 :2 m G11 (q )2 + 2G12 q q + G22 (q )2 = F U + N: (2 40) Dac¼ a în (2 39) înmul¸tim cu g ij ¸si sum¼ am dup¼ a i ob¸tinem ecua¸tia echivalent¼ a :: j m q + j :1 2 11 (q ) +2 j :1 :2 12 q q + j :2 2 22 (q ) = Q1 g 1j +Q2 g 2j ; j = 1; 2: (2 41) Mi¸ scarea pe o suprafa¸ ta tie sub in‡uen¸ ta for¸ tei de greutate ¼ de rota¸ Fie în raport cu reperul Oxyz suprafa¸ta de rota¸tie dat¼ a de ecua¸tiile parametrice x = f (z) cos ; y = f (z) sin ; z = z: Cu nota¸tiile folosite mai sus facem identi…c¼ arile q 1 = ; q 2 = z ¸si avem r0 = f (z) sin i + f (z) cos j r0z = f 0 (z) cos i + f 0 (z) sin j + k; ¼ PUNCTUL MATERIAL SUPUS LA LEGATURI 31 r00 = f (z) cos i f (z) sin j; r00z = f 0 (z) sin i + f 0 (z) cos j; r00zz = f 00 (z) cos i + f 00 (z) sin j; g11 = r0 r0 = f (z)2 ; g12 = r0 r0z = 0; g22 = r0z r0z = 1 + f 0 (z)2 ; 1;11 = r0 r00 = 0; 2;11 1;12 = r0 r00z = f (z) f 0 (z) ; 1;22 = r0 r00zz = 0; = r0z r00 = f (z) f 0 (z) ; 2;21 = r0z r00z = 0; 0 00 0 00 2;22 = r z r zz = f (z) f (z) ; U= r0 r0 r0z r0z = cos i + sin j f 0 (z) k q ; 1 + f 0 (z)2 f (z) f 00 (z) G11 = U r00 = q ; G12 = U r00zz = 0; G22 = q : 1 + f 0 (z)2 1 + f 0 (z)2 S¼ a consider¼ am mai întâi cazul în care axa de simetrie a suprafe¸tei de rota¸tie este vertical¼ a, deci for¸ta de greutate este F = mg = mgk: Ecua¸tiile (2 41) devin :: :: f (z)2 + 2f (z) f 0 (z) z = :: 1 + f 0 (z)2 z f (z) f 0 (z) :2 F r0 = 0; m :2 + f 0 (z) f 00 (z) z = F r0z = m (2 42) g: (2 43) Dup¼ a ce solu¸tiile (t) ¸si z (t) sunt a‡ate se calculeaz¼ a for¸ta de leg¼ atur¼ a (reac¸tiune) din ecua¸tia 1 0 00 (z) :2 f (z) f mgf 0 (z) :2 +q z A= q + N: (2 44) m @q 1 + f 0 (z)2 1 + f 0 (z)2 1 + f 0 (z)2 Se veri…c¼ a imediat c¼ a ecua¸tia (2 42) are integrala prim¼ a (care de altfel poate … ob¸tinut¼ a ¸si din teorema ariilor) : f (z)2 = C: (2 45) 32 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE Înlocuind în (2 43) rezult¼ a :: 1 + f 0 (z)2 z Cf 0 (z) F r0z :2 0 00 + f (z) f (z) z = = m f (z)3 g (2 46) Ecua¸tia (2 46) are integrala prim¼ a :2 1 + f 0 (z)2 z + C2 + 2gz = 2gh; f (z)2 (2 47) care este de fapt integrala prim¼ a a energiei În cazul cilindrului circular drept ca axa vertical¼ a avem f (z) = l: Din : :: (2 45) rezult¼ a = C iar din (2 46) deducem z = g: S¼ a consider¼ am ¸si cazul cilindrului circular drept cu axa de simetrie orizontal¼ a În acest caz putem lua F = mgi iar ecua¸tiile (2 42) ¸si (2 43) vor … înlocuite cu :: + g sin = 0; l :: z = 0: (2 48) (2 49) Ecua¸tia (2 49) are solu¸tia z (t) = v0z t + z0 : Alegând condi¸tiile ini¸tiale v0z = 0; z0 = 0; mi¸scarea are loc în planul Oxy: În acest caz, modelul mecanic prezentat se nume¸ste pendulul matematic ¸si îl vom studia în detaliu într-o subsec¸tiune urm¼ atoare Mi¸ scarea pe o curb¼ a În cazul mi¸sc¼ arii pe curba de ecua¸tii '1 (x; y; z; t) = 0; '2 (x; y; z; t) = 0 din legea a doua a lui Newton ob¸tinem ma = F + N unde F este rezultanta for¸telor active iar for¸ta de leg¼ atur¼ a (reac¸tiune) este N= 1 grad'1 + 2 grad'2 : ¼ PUNCTUL MATERIAL SUPUS LA LEGATURI 33 În cazul unei curbe …xe, considerând triedrul lui Frenet ata¸sat curbei ; ; ; notând F =F +F +F ; N =N +N avem : mv = F ; mv 2 (2 50) =F +N ; 0=F +N : (2 51) (2 52) Dac¼ a ecua¸tia curbei este dat¼ a sub form¼ a parametric¼ a r = r (s) ; s reprezentând abscisa curbilinie, deducem :: ms = F : s; s; t : : Integrând aceast¼ a ecua¸tie diferen¸tial¼ a cu condi¸tiile ini¸tiale s (t0 ) = s0 ; s (t0 ) = v0 se ob¸tine legea de mi¸scare pe curba dat¼ a s = s (t) : În cazul în care F depinde doar de pozi¸tia punctului, avem din teorema energiei cinetice (2 19) : Z t Z t Z s 1 1 : 2 2 mv mv = F vdt = F sdt = F (s) ds; (2 53) 2 2 0 t0 t0 s0 : de unde deducem v = v (s) : Având pe v = s, putem a‡a timpul printr-o cuadratur¼ a Z s ds t t0 = ; s0 v (s) de unde ob¸tinem apoi s = s (t) : Dac¼ a func¸tia F are poten¸tialul U; atunci F v = grad U dr dU = dt dt ¸si din (2 53) deducem integrala prim¼ a a energiei T U = const: 34 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE Pendulul matematic Consider¼ am un reper cartezian Oxyz cu axa Ox având direc¸tia ¸si sensul verticalei descendente Ne propunem s¼ a studiem mi¸scare unui punct material greu M pe cercul vertical x2 + y 2 = l2 ; z = 0: Pentru simplitate vom omite s¼ a mai scriem a treia variabil¼ a spa¸tial¼ a Normala la cerc este de forma N = xi + yj iar ecua¸tiile de mi¸scare sunt :: :: mx = x + mg; my = y: \ Fie = (Ox; OM ) Trecând la coordonate polare x = l cos ; y = l sin (…gura 2 2) ecua¸tiile de mi¸scare devin :: ml :2 sin :: ml cos :2 cos = l cos + mg; (2 54) sin = l sin : Figure 2 2: Pendulul matematic Prin eliminarea lui se ob¸tine ecua¸tia pendulului matematic :: + g sin = 0: l :: Eliminând g¼ asim valoarea lui = ml :2 mg cos ; l (2 55) ¼ PUNCTUL MATERIAL SUPUS LA LEGATURI 35 ¸si de aici m¼ arimea algebric¼ a a reac¸tiunii :2 N = mg cos + ml : (2 56) Putem reg¼ asi mai u¸sor rela¸tiile (2 55), (2 56) scriind ecua¸tiile de mi¸scare în proiec¸tie pe axele triedrului lui Frenet Tinând ¸ cont c¼ a tangenta ¸si normala la curb¼ a sunt = sin i + cos j; iar viteza : = cos i + sin j; : : v = xi + yj = l sin : i + l cos j; ecua¸tiile (2 50), (2 51) devin :: ml = ml :2 = mg sin ; mg cos + N; adic¼ a tocmai ecua¸tiile (2 55), (2 56) Ecua¸tia pendulului matematic poate … ob¸tinut¼ a ¸si din teorema momentului cinetic Momentul cinetic fa¸ta¼ de O este KO = r mv = ml2 cos i + sin j : sin : i + cos : j = ml2 k iar momentul rezultant al for¸telor de greutate ¸si reac¸tiune MO = r N +r mg = r r + l cos i + sin j mgi = mgl sin k: Aplic¼ am teorema momentului cinetic ¸si g¼ asim ecua¸tia pendulului matematic: : :: KO = M O ) + g sin = 0: l : Înmul¸tind (2 55) cu l2 ¸si integrând g¼ asim integrala prim¼ a l2 :2 2gl cos = h = const: (2 57) :2 1 1 Cum energia cinetic¼ a a punctului material este T = mv 2 = ml2 iar 2 2 poten¸tialul for¸tei de greutate este U = mgx = mgl cos ; integrala prim¼ a (2 57) este chiar integrala prim¼ a a energiei Tipuri de mi¸ sc¼ari ale punctului material greu pe cerc 36 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE Mi¸scarea punctului material greu pe un cerc vertical va … determinat¼ a de condi¸tiile ini¸tiale : : (0) = 0 ; (0) = 0 : (2 58) M¼ arimea vitezei ini¸tiale ¸si abscisa pozi¸tiei ini¸tiale sunt : v0 = l 0 ; x0 = l cos 0: Integrala prim¼ a a energiei (2 57) se expliciteaz¼ a l2 :2 2gl cos = v02 2gx0 : (2 59) Fie xA abscisa de…nit¼ a prin xA = x0 v02 ; 2g astfel încât integrala prim¼ a se scrie sub forma l2 :2 = 2gl(cos xA ): (2 60) Aceast¼ a rela¸tie ne arat¼ a c¼ a în cursul mi¸sc¼ arii xA x l: 1 Cazul l 0: Din ecua¸tia de mai sus, în care se ia semnul " + "; rezult¼ a c¼ a : va cre¸ste pân¼ a la valoarea :când se anuleaz¼ a Luând apoi semnul " " în fa¸ta radicalului, se vede c¼ a va descre¸ste pân¼ a la valoarea , când iar¼ a¸si ¼ PUNCTUL MATERIAL SUPUS LA LEGATURI : 37 : se anuleaz¼ a etc Intervalul de timp în care ¸si revin la aceea¸si valoare este perioada de oscila¸tie s Z d l r = : (2 63) g 2 2 sin sin 2 2 Notând sin = u sin ; 2 2 sub forma integralei eliptice s Z l 1 du p =4 2 g 0 (1 u ) (1 se exprim¼ a k 2 u2 ) ; (2 64) (forma normal¼ a Legendre), unde parametrul k 2 = sin2 depinde de condi¸tiile ini¸tiale Pentru calculul numeric al lui p (1 1 k 2 u2 ) =1+ 2 2g (x0 + l) ¸si punctul descrie complet cercul în acela¸si : sens, impus de semnul lui 0 Ecua¸tia diferen¸tial¼ a care descrie mi¸scarea este : l = s 2g (l xA ) 1 unde k12 sin2 2 2l 0 ¸si q (r) > 0 pentru r R: Cum la : momentul ini¸tial r = R, deducem în acest caz c¼ a r nu se anuleaz¼ a, r¼ amâne tot timpul pozitiv ¸si distan¸ta de la centrul P¼ amântului la proiectil cre¸ste în timp Dac¼ a v0 0 ¸si limr!1 q (r) R; r se anuleaz¼ a într-un punct (punctul cel mai înalt al traiectoriei) unde avem : r = rmax ; r = 0 ¸si din (4 13) deducem ecua¸tia algebric¼ a 2gR 2 v02 rmax 2gR2 rmax + R2 v02 cos2 = 0: 58 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE din care rezult¼ a gR2 rmax = 2gR p 2 : v02 Cu condi¸tia rmax R din solu¸tiile de mai sus se re¸tine cea cu semnul ’+’ Rezultatele ob¸tinute aici vor … reg¼ asite în sec¸tiunea care urmeaz¼ a mar = F 4 1 3 mat maC : (4 14) Mi¸ scarea corpurilor grele For¸ ta de greutate Vom considera în cotinuare sistemul geocentric, adic¼ a un reper (presupus iner¸tial) având originea O1 în centrul P¼ amântului, axa O1 z1 orientat¼ a dup¼ a linia polilor (adic¼ a linia ce une¸ste polul N cu polul S) iar axele O1 x1 ¸si O1 y1 în plan ecuatorial orientate spre dou¼ a stele presupuse …xe Rota¸tia : P¼ amântului este ! = !k 1 cu k 1 versorul axei O1 z1 ¸si practic ! = 0 ) ! = const: S¼ a alegem un reper mobil cartezian solidar legat de P¼ amânt astfel încât axa O1 z s¼ a coincid¼ a cu axa O1 z1 iar axele O1 x ¸si O1 y s¼ a …e în planul ecuatorial (…gura 4 1, a)) Fie la suprafa¸ta P¼ amântului punctul material M de mas¼ a m: Asupra punctului ac¸tioneaz¼ a din partea P¼ amântului for¸ta de atrac¸tie newtonian¼ a (al c¼ arei suport trece prin centrul P¼ amântului) maT a mai trebuie ad¼ augat un termen pentru a ob¸tine greutatea) (nu mg pentru c¼ 0 precum ¸si for¸te de alt¼ a natur¼ a a c¼ aror rezultant¼ a o not¼ am cu F : Rezult¼ a deci 0 (4 15) mar = F + maT ! (! r) 2m! v r : ¸tinând cont c¼ a ! (! r) = !2P M = ! 2 d; unde P este proiec¸tia lui M pe O1 X3; avem 0 mar = F + mg 2m! vr ; (4 16) cu g = aT + ! 2 d: (4 17) Atrac¸tia P¼ amântului maT plus for¸ta centrifug¼a m! 2 d dau for¸ta de greutate mg unde g reprezint¼ a accelera¸tia gravita¸tional¼a Suportul lui g nu trece neap¼ arat prin centrul O1 dar este în planul meridian ce trece prin punctul M: ¼ MECANICA ¼ TERESTRA ¼ MI¸ SCAREA RELATIVA 59 Figure 4 1: Devierea …rului cu plumb 4 1 4 Devierea …rului cu plumb Varia¸ tia accelera¸ tiei gravita¸ tionale la suprafa¸ ta P¼ amântului S¼ a studiem echilibrul …rului cu plumb (…gura 4 1, b)) Asupra bilei metalice 0 în echilibru ac¸tioneaz¼ a for¸ta de tensiune în …r F = T ¸si greutatea mg: Condi¸tia de echilibru ar = 0; v r = 0 conduce la rela¸tia T + mg = 0: Vectorul g se g¼ ase¸ste în planul meridian al punctului M: Consider¼ am c¼ a aT este riguros dirijat spre centrul O1 al P¼ amântului ¸si not¼ am cu unghiul dintre raza O1 M , suport al lui aT ¸si verticala locului (suportul lui g) m¼ asoar¼ a devierea …rului cu plumb Not¼ am cu latitudinea punctului M (unghiul dintre O1 M ¸si planul ecuatorial Vectorul T de sens opus lui g face cu planul ecuatorial unghiul + : Proiectând rela¸tia (4 17) pe verticala ¸si orizontala locului avem g = aT cos 0 = aT sin ! 2 d cos ( + ) ; ! 2 d sin ( + ) ; (4 18) cu d = R cos ; unde R este raza P¼ amântului Cum sin (4 19) este mic, cu aproxima¸tiile ; sin ( + ) sin ; 60 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE din (4 18) ¸si (4 19) deducem c¼ a = ! 2 R sin 2 : 2aT A¸sadar devierea maxim¼ a a …rului cu plumb de la direc¸tia perpenducular¼ a pe suprafa¸ta P¼ amântului se ob¸tine pentru = : La aceast¼ a latitudine avem 4 6 grade sexagesimale 4 1 5 Ecua¸ tiile de mi¸ scare la suprafa¸ ta P¼ amântului Vom considera în cele ce urmeaz¼ a c¼ a g este practic constant ¸si dirijat spre centrul P¼ amântului Consider¼ am la suprafa¸ta P¼ amântului reperul Oxyz solidar legat de P¼ amânt ales astfel (…gura 4 2, a)): Oxeste dirijat pe tangenta la paralela ce trece prin punctul O cu sensul pozitiv spre est, Oy este dirijat pe tangenta la meridianul ce trece prin O cu sensul pozitiv spre nord iar Oz este verticala ascendent¼ a Evident avem ! = ! cos i2 + ! sin i3 : iar în cea sudic¼ a 2 2 Proiec¸tiile ecua¸tiei (4 14) pe axele reperului Oxyz sunt În emisfera nordic¼ a0 :: : : mx = F10 2m! z cos y sin ; :: : my = F20 2m! x sin ; :: : mz = F30 mg + 2m! x cos : 4 1 6 0: (4 20) Pendulul lui Foucault O prob¼ a mecanic¼ a a rota¸tiei P¼ amântului o aduce celebra experien¸ta¼ a lui L Foucault efectuat¼ a în 1851 la Observatorul Astronomic din Paris, ¸si repetat¼ a în mod public în 1852 la Pantheonul din Paris Fie un pendul format dintr-o bil¼ a metalic¼ a M ata¸sat¼ a de un …r inextensibil al c¼ arui cap¼ at este …xat în O În acest caz r 0 F = T = N ; r = OM ; l = OM ; l iar ecua¸tiile (4 20) devin :: : : N xl1 2m! z cos y sin :: : my = N yl 2m! x sin ; :: : mz = N zl mg + 2m! x cos ; mx = ; (4 21) ¼ MECANICA ¼ TERESTRA ¼ MI¸ SCAREA RELATIVA 61 Figure 4 2: Repere neiner¸tiale la suprafa¸ta Pamântului cu x2 + y 2 + z 2 = l 2 : : : (4 22) : Înmul¸tind ecua¸tiile (4 21) respectiv cu x; y; z; adunând, integrând în raport cu t ¸si ¸tinând cont de (4 22) ob¸tinem integrala prim¼ a a energiei :2 :2 :2 x + y + z = 2g (z + h) : Din (4 22) deducem r z= l 1 x2 + y = l2 l 1 (4 23) x2 + y 2 + ::: : 2l2 x y ; vor … mici, de l l ordinul lui 0 0: De…nim o conic¼a drept mul¸timea punctelor M din plan pentru care kOM k = e = const: kM M 0 k (8 1) unde M 0 reprezint¼ a proiec¸tia lui M pe dreapta (d) iar e > 0 reprezint¼ a excentricitatea conicei Punctul O se nume¸ste focarul conicei iar dreapta (d) este directoarea conicei Dac¼ a e 1 este o hiperbol¼a Notând cu (x; y) coordonatele carteziene ale punctului M , (8 1) revine la x2 + y 2 = e2 (x , r2 1 k)2 , r2 = e2 (r cos e2 cos 2 de unde deducem r1;2 = + 2e2 kr cos ek (e cos e2 cos 2 k)2 , e2 k 2 = 0; 1) ; 1 sau echivalent p p ; r2 = ; (8 2) 1 + e cos 1 e cos unde p = ke este parametrul conicei a) Cazul elipsei Tinând ¸ cont c¼ a e 0 deducem c¼ a din cele dou¼ a solu¸tii (8 2) putem alege numai solu¸tia r1 = r= p : 1 + e cos 105 (8 3) 106 Notând cu a semiaxa mare a elipsei 2a = p 1 1 + 1+e 1 e = 1 2p : e2 Distan¸ta dintre focarele O ¸si O0 ( 2c; 0) este 2c = r ( ) r (0) = 2pe : 1 e2 (8 4) Din (8 3) ¸si (8 4) deducem e= c ; p=a 1 a a2 e2 = c2 a = b2 ; a unde b este semiaxa mic¼ a a elipsei b) Cazul parabolei Cum e = 1, avem ecua¸tia parabolei în coordonate polare q p = ; (8 5) r= 1 + cos cos2 2 unde p = 2q = r(0): (8 6) c) Cazul hiperbolei (e > 1) Din (8 2) rezult¼ a c¼ a hiperbola are dou¼ a ramuri r1 ( ) = p ; 1 + e cos r2 ( ) = 1 2 p ; e cos 1 e arccos 2 arccos 1 e 1 e ; arccos ; arccos 1 e ; : Distan¸ta între punctele în care cele dou¼ a ramuri ale hiperbolei intersecteaz¼ a axa Ox este 2a = r2 (0) r1 (0) = p e p 2p = 2 : e+1 e 1 1 Distan¸ta dintre focarele hiperbolei este 2c = r2 (0) + r1 (0) = de unde deducem e= p e 1 + p 2pe = 2 ; e+1 e 1 c b2 ; p = ; b2 = c2 a a a2 : 107 8 0 6 Deducerea legii de atrac¸ tie newtonian¼ a din legile lui Kepler În a doua jum¼ atate a secolului al XVI-lea, astronomul danez Tycho Brahe reu¸se¸ste s¼ a strâng¼ a o mare cantitate de observa¸tii foarte precise referitoare la pozi¸tiile planetei Marte În anii 1609 ¸si 1619 Johannes Kepler, prelucrînd observa¸tiile astronomice ale lui Tycho Brahe a enun¸tat în urma urm¼ atoarele legi: 1 Planetele asimilate cu puncte materiale descriu în jurul Soarelui traiectorii plane s¸i anume elipse, Soarele …ind plasat în focar 2 Mi¸ scarea are loc respectându-se legea (teorema) ariilor 3 Raportul dintre a3 ; cubul semiaxei mari a elipsei s¸i T 2 ; p¼atratul timpului de revolu¸tie este acela¸ si pentru toate planetele Primele dou¼ a legi au fost enun¸tate în lucrarea Astronomia nova (1609) iar a treia a fost dat¼ a în lucrarea De Harmonices mundi (1619) Consider¼ am mi¸scarea în planul Oxy; Soarele …ind plasat în O: Vectorul : : : de pozi¸tie al unei planete …ind r = ri ; viteza sa va … v = ri + r i : : Deducem de aici, în virtutea legii ariilor, c¼ a r v = r2 i3 = Ci3 : Fie F 108 for¸ta exercitat¼ a de Soare asupra planetei Teorema momentului cinetic ne furnizeaz¼ a rela¸tiile : i d h d Ci3 = 0 ) r k F r F = m r r =m dt dt deci F este o for¸ta¼ central¼ a F = Fi : Fie a ¸si b semiaxele elipsei ce reprezint¼ a, în virtutea primei legi a lui Kepler, traiectoria planetei ¸si p a2 b2 b2 e= 0 ) e > 1 deci traiectoria este o hiperbol¼ a r0 Constanta 2 ; r0 h = v02 este constanta care apare aplicând legea conserv¼ arii energiei la problema studiat¼ a: T 8 2 1 m 2 U= 2 r v2 = h = const: Traiectorii (orbite) eliptice Sateli¸ tii arti…ciali ai p¼ amântului S¼ a presupunem c¼ a datele ini¸tiale sunt astfel încât h 1 0e-4); E0=E; E=E-(E-e*sin(E)-Mt)/(1-e*cos(E)); end f=E; end PROBLEMA LUI NEWTON 119 Figure 8 3: Mi¸scarea pe orbita eliptic¼ a Utilizând formulele (8 26), am conceput programul de anima¸tie MATLAB miscare pe elipsa m care prezint¼ a mi¸scarea pe orbita eliptic¼ a în jurul unui pol atractiv %*************miscare pe elipsa m********* t=0:T/36:T; s=size(t); n=2*pi/T; for i=1:s(2), u(i)=kepler(e,n*t(i)); csi(i)=a*(cos(u(i))-e); eta(i)=a*sqrt(1-e^2)*sin(u(i)); for j=1:N, cs(j,i)=j*csi(i)/N; et(j,i)=j*eta(i)/N; end end M0=moviein(s(2)); for i=1:s(2), plot(csi,eta,’LineWidth’,2); hold on; plot(cs(:,i),et(:,i)); 120 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE plot(csi(i),eta(i),’o’,’LineWidth’,2) M0(:,i)=getframe; end movie(M0); În …guea 8 3 prezent¼ am mi¸scarea pe elips¼ a pentru valorile T = 1; a = 1; e = 0:4; = 0: 8 3 2 Traiectoriile în cazul mi¸ sc¼ arii absolute Consider¼ am mi¸scarea corpurilor P1 (având masa M ¸si vectorul de pozi¸tie r1 ) ¸si P2 (având masa M ¸si vectorul de pozi¸tie r2 ) relativ la un reper iner¸tial cu originea în centrul de mas¼ a Avem a¸sadar sistemul algebric liniar M r1 + mr2 = 0; r1 r2 = r; cu solu¸tia r1 = 1 M=m r; r2 = r; 1 + M=m 1 + M=m (8 32) r …ind calculat mai sus În programul doua corpuri m scris în MATLAB calcul¼ am r1 ¸si r2 : %***********doua corpuri********* t=0:T/36:T; s=size(t); n=2*pi/T; for i=1:s(2), u(i)=kepler(e,n*t(i)); csi(i)=a*(cos(u(i))-e); eta(i)=a*sqrt(1-e^2)*sin(u(i)); csi1(i)=-csi(i)*m/(M+m); csi2(i)=csi(i)*M/(M+m); eta1(i)=-eta(i)*m/(M+m); eta2(i)=eta(i)*M/(M+m); for j=1:N, cs1(j,i)=j*csi1(i)/N; et1(j,i)=j*eta1(i)/N; cs2(j,i)=j*csi2(i)/N; et2(j,i)=j*eta2(i)/N; end end PROBLEMA LUI NEWTON 121 Figure 8 4: Traiectoriile celor dou¼ a corpuri M0=moviein(s(2)); for i=1:s(2), plot(csi1,eta1,’LineWidth’,2); hold on; plot(csi2,eta2,’LineWidth’,2); if (rem(i,4)==0) plot(cs1(:,i),et1(:,i),’LineWidth’,0 5); plot(cs2(:,i),et2(:,i),’LineWidth’,2); end plot(csi1(i),eta1(i),’o’,’LineWidth’,2) plot(csi2(i),eta2(i),’o’,’LineWidth’,2) M0(:,i)=getframe; end movie(M0); În …gura 8 4 sunt prezentate traiectoriile corpurilor P1 ¸si P2 pentru diverse valori ale raportului M=m S-au considerat valorile valorile T = 1; a = 1; e = 0:4; = 0: 122 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE Figure 8 5: Planul ecliptic Ilustra¸tie a lui Tim Gunthe, National Geographic 8 3 3 Sistemul Solar (material preluat din Wikipedia) Sistemul solar este format din Soare împreun¼ a cu sistemul s¼ au planetar (care cuprinde opt planete împreun¼ a cu sateli¸tii lor naturali) ¸si alte obiecte non-stelare Sistemul este situat într-unul dintre bra¸tele exterioare ale galaxiei Calea Lactee (mai precis în Bra¸tul Orion), galaxie care are cca 200 de miliarde de stele El s-a format acum 4,6 miliarde de ani, ca urmare a colapsului gravita¸tional al unui gigant nor molecular Cel mai masiv obiect este steaua central¼ a - Soarele, al doilea obiect ca mas¼ a …ind planeta Jupiter Cele patru planete interioare mici, Mercur, Venus, P¼ amântul ¸si Marte, numite ¸si planete terestre / planete telurice, sunt compuse în principal din roci ¸si metal Cele patru planete exterioare, numite ¸si gigan¸ti gazo¸si, sunt mult mai masive decât cele terestre Cele mai mari dou¼ a planete, Jupiter ¸si Saturn, sunt compuse în principal din hidrogen ¸si heliu; cele dou¼ a planete mai îndep¼ artate, Uranus ¸si Neptun, sunt compuse în mare parte din substan¸te cu o temperatur¼ a de topire relativ ridicat¼ a (comparativ cu hidrogenul ¸si heliu), numite ghe¸turi, cum ar … apa, amoniacul ¸si metanul Ele sunt denumite „gigan¸ti de ghea¸ta¼” (termen distinct de cel de „gigant gazos”) Toate planetele au orbite aproape circulare dispuse într-un disc aproape plat numit plan ecliptic Sistemul solar prezint¼ a câteva regiuni unde se a‡a¼ diferite obiecte mici Centura de asteroizi, situat¼ a între Marte ¸si Jupiter, este similar¼ a din punct de vedere al compozi¸tiei cu planetele terestre, deoarece o mare parte dintre obiecte sunt compuse din roc¼ a ¸si metal Dincolo de orbita lui Neptun se a‡a¼ centura Kuiper ¸si discul împr¼ a¸stiat; multe dintre obiectele transneptuniene sunt în mare parte compuse din ghe¸turi Printre aceste obiecte, de la câteva zeci pân¼ a la mai mult de zece mii PROBLEMA LUI NEWTON 123 de obiecte pot … su…cient de mari pentru a … fost rotunjite de propria gravita¸tie Astfel de obiecte sunt denumite planete pitice Planetele pitice identi…cate pân¼ a în prezent includ asteroidul Ceres ¸si obiectele transneptuniene: Pluto, Eris, Haumea ¸si Makemake În plus, în aceste dou¼ a regiuni se a‡a¼ diferite alte corpuri mici, cum ar … comete, centauri ¸si materie interplanetar¼ a S ¸ase planete, cel pu¸tin trei planete pitice ¸si multe alte corpuri mici au sateli¸ti naturali care se rotesc în jurul lor Fiecare planet¼ a exterioar¼ a este înconjurat¼ a de inele planetare alc¼ atuite din praf ¸si alte obiecte mici Vântul solar (un ‡ux de plasm¼ a de la Soare) creeaz¼ a în mediul interstelar o bul¼ a cunoscut¼ a ca heliosfer¼ a, care se extinde pân¼ a la marginea discului împr¼ a¸stiat La limita sa exterioar¼ a se a‡a¼ Norul lui Oort, care reprezint¼ a doar un câmp de resturi r¼ amase dup¼ a crearea planetelor, …ind considerat a … sursa pentru cometele cu perioad¼ a lung¼ a El se întinde pân¼ a la o distan¸ta¼ aproximativ de o mie de ori mai mare decât heliosfera Heliopauza este punctul în care presiunea vântului solar este egal¼ a cu presiunea opus¼ a a vântului interstelar 8 4 Formulele mi¸ sc¼ arii parabolice Cometele pot avea în raport cu Soarele orbite parabolice, eliptice sau hiperbolice Cele cu orbite eliptice se numesc comete periodice O astfel de comet¼ a este cometa Halley care a trecut la periheliu ultima dat¼ a în 1986, perioada ei …ind de aproximativ 75 de ani În cazul mi¸sc¼ arii parabolice (…gura 8 6) avem e = 1 ¸si ecua¸tia orbitei parabolice în coordonate polare r= p = 1 + cos # unde q= q cos2 # 2 ; (8 33) p = kP1 k 2 (8 34) reprezint¼ a distan¸ta perihelic¼ a (pericentric¼ a) Din teorema ariilor scris¼ a sub forma d# r2 =C dt ¸si din rela¸tia (8 33) deducem Cdt = r2 d# = q2 cos4 # 2 d# = 2q 2 1 + tan2 # 2 d tan # 2 ; (8 35) 124 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE Figure 8 6: Elementele orbitei parabolice PROBLEMA LUI NEWTON de unde, integrând de la t = C (t )= Z tan # 2 125 la momentul curent t ob¸tinem 2q 2 1 + tan2 0 # 2 d tan Cum p= C 2 # 2 = 2q 2 tan # 1 # : + tan3 2 3 2 (8 36) ; (vezi rela¸tia (8 24)), din (8 36) ob¸tinem # 1 # tan + tan3 = 2 3 2 r 2q 3 (t ): (8 37) Rezolv¼ and ecua¸tia algebric¼ a (8 37) ¸si folosind (8 36) putem calcula coordonatele ( ; ) ale corpului P2 : = r cos # = r cos2 # 2 sin2 # 2 =q 1 = r sin # = 2r cos 8 4 1 = r cos2 tan2 # 2 # 2 1 ; # # # sin = 2q tan : 2 2 2 tan2 # 2 = (8 38) (8 39) Metoda bisec¸ tiei Metoda bisec¸tiei const¼ a în reducerea intervalului de c¼ autare prin înjum¼ at¼ a¸tirea repetat¼ a ¸si selectarea subintervalului în care se g¼ ase¸ste r¼ ad¼ acina De aceea algoritmul este logaritmic n(O(nlog{n})n) Algoritmul este similar cu c¼ autarea binar¼ a Ca metoda s¼ a func¸tioneze trebuie ca func¸tia, c¼ areia vrem s¼ a-i g¼ asim r¼ ad¼ acina, s¼ a aib¼ a proprietatea lui Darboux (de exemplu s¼ a …e continu¼ a) Dac¼ a f (o func¸tie oarecare) este continu¼ a pe intervalul [a; b] ¸si f (a) ¸si f (b) au semne opuse, adic¼ a f (a) f (b) 1 0e-4) c=(a+b)/2; if ((a+a^3/3-tl)*(c+c^3/3-tl)>0) a=c; else b=c; end end f=a; 8 4 2 Mi¸ scarea parabolic¼ a relativ¼ a Vom folosi coordonate adimensionale Adimensionaliz¼ am timpul prin rela¸tia r t = (t ); 2q 3 ¸si coordonatele spa¸tiale prin rela¸tiile = ; q = q : Pentru a studia mi¸scarea corpului P2 în jurul polului (corpului) atractiv P1 în intervalul adimensional t 2 [ T; T ] ; rezolv¼ am ecua¸tia polinomial¼ a de T3 T3 gradul 3(8 37), folosind metoda bisec¸tiei în intervalul T ;T + : 3 3 Cu formulele (8 38), (8 39) calcul¼ am (t ); (t ): Cu programul MATLAB miscare pe parabola m, prezent¼ am (…gura8 7) în coordonate adimensionale pozi¸tiile pe orbita parabolic¼ a (¸si razele vectoare corespunz¼ atoare) pe care corpul P2 le are la momente echidistante %*********miscare pe parabola m******** T=20;N=40; t=-T:T/101:T; s=size(t); for i=1:s(2), PROBLEMA LUI NEWTON 127 a=-T-T^3/3; b=-a; sol(i)=bisectie(a,b,t(i)); csi(i)=1-sol(i)^2; eta(i)=2*sol(i); for j=1:N, cs(j,i)=j*csi(i)/N; et(j,i)=j*eta(i)/N; end end M=moviein(s(2)); for i=1:s(2), plot(csi,eta,’LineWidth’,2); hold on; if (rem(i,6)==0) plot(cs(:,i),et(:,i)); plot(csi(i),eta(i),’o’,’LineWidth’,2) end M(:,i)=getframe; end movie(M); 8 4 3 Mi¸ scarea parabolic¼ a absolut¼ a Folosim formulele (??) pentru a studia mi¸scarea absolut¼ a în raport cu un reper iner¸tial cu originea în centrul de mas¼ a al sistemului format din cele dou¼ a corpuri Cu programul MATLAB de mai jos calcul¼ am orbitele corpurilor (punctelor) P1 ¸si P2 în mi¸scarea absolut¼ a ¸si prezent¼ am (…gura 8 8) pozi¸tiile corpurilor ¸si razele vectoare corespunz¼ atoare la momemnte echidistante Considerând diverse valori ale raportului dintre mase M=m: T=20;N=40; t=-T:T/101:T; s=size(t); for i=1:s(2), a=-T-T^3/3; b=-a; sol(i)=bisectie(a,b,t(i)); csi(i)=1-sol(i)^2; eta(i)=2*sol(i); csi1(i)=-csi(i)*m/(M+m); 128 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE Figure 8 7: Mi¸scarea relativ¼ a pe orbita parabolic¼ a csi2(i)=csi(i)*M/(M+m); eta1(i)=-eta(i)*m/(M+m); eta2(i)=eta(i)*M/(M+m); for j=1:N, cs1(j,i)=j*csi1(i)/N; et1(j,i)=j*eta1(i)/N; cs2(j,i)=j*csi2(i)/N; et2(j,i)=j*eta2(i)/N; end end M0=moviein(s(2)); for i=1:s(2), plot(csi1,eta1,’LineWidth’,2); hold on; plot(csi2,eta2,’LineWidth’,2); if (rem(i,8)==0) plot(cs1(:,i),et1(:,i)); plot(cs2(:,i),et2(:,i)); plot(csi1(i),eta1(i),’o’,’LineWidth’,2) plot(csi2(i),eta2(i),’o’,’LineWidth’,2) PROBLEMA LUI NEWTON 129 Figure 8 8: Mi¸scarea absolut¼ a în cazul orbitelor parabolice end M0(:,i)=getframe; end movie(M0); 8 5 8 5 1 Pozi¸ tia orbitei în spa¸ tiu Rela¸ tia matriceal¼ a de transformarea sistemelor de coordonate Fie P1 planul orbitei corpului P2 ; atras de P1 cu for¸ta de atrac¸tie universal¼ a Consider¼ am ¸si axa P1 ? P1 astfel încât P1 s¼ a …e un reper ortogonal Drept ax¼ a P1 lu¼ am axa P1 care trece prin pericentrul P1 se mai nume¸ste ¸si sistemul de coordonate orbitale carteziene Al¼ aturi de acest reper vom considera un reper ortogonal de referin¸t¼a P1 xyz Fie dreapta N N 0 (linia nodurilor), intersec¸tia planului P1 al orbitei cu planul fundamental P1 xy (…gura 8 9) Sensul pozitiv al liniei nodurilor este cel al semidreptei P1 N pe care punctul P2 o intersecteaz¼ a trecând din domeniul z 0 130 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE Figure 8 9: Pozitia planului orbitei PROBLEMA LUI NEWTON 131 Unghiul = (P2\ x; P2 N ) 2 00 ; 3600 poart¼ a numele de longitudinea nodului ascendent ¸si se m¼ asoar¼ a de la axa P1 x spre axa P1 y Fie unghiul 0 0 \ i = (P1 ; P1 z) 2 0 ; 180 numit înclinarea orbitei (sau în¼al¸timea orbitei ) Dac¼ a 00 < i < 900 , mi¸scarea este direct¼a iar dac¼ a 900 < i < 1800 , mi¸scarea este retrograd¼a : Vom trece de la reperul P1 xyz la un reper ortogonal P1 x1 y1 z1 al c¼ arui plan fundamental este chiar planul orbitei (deci P1 = P1 z1 ) ¸si a c¼ arui ax¼ a P1 x1 coincide cu semiaxa nodurilor P1 N În planul P1 x1 y1 deplasarea punctului P2 are loc în sensul P2 x1 ! P2 y1 Transformarea de coordonate (x1 ; y1 ; z1 ) ! (x; y; z) va avea loc în dou¼ a etape În prima etap¼ a se trece de la reperul P1 xyz la reperul ajut¼ ator P1 x1 y 0 z, ob¸tinut din P1 xyz prin rotirea în jurul axei P1 z cu unghiul Formula de trecere este 10 1 0 1 0 x1 cos sin 0 x @ y A = @ sin cos 0 A @ y0 A : (8 40) 0 0 1 z z În etapa a doua se trece de la reperul P1 x1 y 0 z la P1 x1 y1 z1 prin rota¸tia de unghi i în jurul axei P1 x1 Formula de trecere este 10 0 1 0 1 x1 x1 1 0 0 @ y 0 A = @ 0 cos i sin i A @ y1 A : (8 41) 0 sin i cos i z z1 Prin compunerea transform¼ arilor (8 40) ¸si (8 41) ob¸tinem trecerea de la reperul P1 xyz la reperul P1 x1 y1 z1 : 0 1 0 10 1 x cos cos i sin sin i sin x1 @ y A = @ sin cos i cos sin i cos A @ y1 A : (8 42) z 0 sin i cos i z1 8 5 2 Elementele orbitei Fie r (t) = P1 P2 = (x; y; z) ¸si r0 = (x0 ; y0 ; z0 ) valuarea vectorului de pozi¸tie P1 P2 la un moment t0 Pentru a rezolva problema lui Cauchy pentru sistemul de ecua¸tii diferen¸tiale de ordinul doi, trebuie s¼ a cunoa¸stem ¸si viteza a cum ¸stim din teorema ariilor la momentul ini¸tial t0 ; r0 = x0 ; y 0 ; z 0 : Dup¼ r r = r0 r0 = C = const: 132 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE iar mi¸scarea are loc pe o orbit¼ a plan¼ a, a‡at¼ a în planul P1 x1 y1 : Coordonatele vectorului C în raport cu reperul P1 x1 y1 z1 vor … C = (0; 0; C) : (8 43) Notând cu (Cx ; Cy ; Cz ) coordonatele lui C în reperul P1 xyz, avem în virtutea transform¼ arii (8 42): 1 0 10 1 0 cos cos i sin sin i sin 0 Cx @ Cy A = @ sin cos i cos sin i cos A @ 0 A : (8 44) 0 sin i cos i C Cz Din (8 44), cu teorema ariei ob¸tinem yz z y = C sin i sin ; (8 45) zx xz = (8 46) xy y x = C cos i; C sin i cos ; (8 47) rela¸tii care permit calcularea constantelor C; i; ; cunoscând pozi¸tia ¸si viteza ini¸tial¼ a a punctului P2 În sistemul de coordonate P2 x1 y1 z1 ; mi¸scarea având loc în planul P2 x1 y1 ; z1 = 0; z 1¸si teorma ariilor conduce la x1 y 1 y1 x1 = C; (8 48) expresie care, dup¼ a introducerea coordonatelor polare r = kP1 P2 k ; \ r; P1 x1 prin formulele = x1 = r cos ; (8 49) y1 = r sin ; (8 50) z1 = 0; (8 51) conduce la r2 = C: Unghiul polar se nume¸ste argumentul de latitudine Din (8 42), (8 49), (8 50) ¸si (8 51) rezult¼ a coordonatele carteziene ale lui P2 în reperul P1 xyz ca func¸tii de coordonatele polare din planul P1 x1 y1 : x = r (cos cos sin sin cos i) ; y = r (cos sin + sin cos cos i) ; z = r sin sin i: PROBLEMA LUI NEWTON 133 Componentele vitezei relative sunt x + r ( sin cos cos sin r y y = r + r ( sin sin + cos cos r z z = r + r cos sin i: r x = r cos i) ; cos i) ; Am ob¸tinut astfel coordonatele s¸i vitezele parametrice ale mi¸ sc¼arii relative În locul sistemului P1 x1 y1 z1 putem folosi sistemul de coordonate orbitale carteziene P1 Fie longitudinea pericentrului (periastrului) ! = (P1\ N; P1 ): Reperul P1 se ob¸tine din P1 x1 y1 z1 prin rotite cu unghiul ! în jurul axei P1 = P1 z1 Avem deci 10 1 0 1 0 x1 cos ! sin ! 0 @ y1 A = @ sin ! cos ! 0 A @ A 0 0 0 1 z1 ¸si apoi cu (8 42) putem exprima coordonatele x; y; z în func¸tie de ¸si : Reamintim c¼ a am notat cu momentul trecerii lui P1 la pericentru (periastru), cu e excentricitatea orbitei ¸si cu q = kP1 k, distan¸ta de la P1 la pericentrul Elementele ( ; !; i; e; q; ) poart¼ a numele de elementele orbitei Cunoscând elementele orbitei ¸si valoarea constantei C a ariilor se poate calcula r (t) : 8 6 Calculul efemeridelor pentru planete, comete ¸ si sateli¸ ti Not¼ am S = P1 ; P = P2 Nota¸tia S o folosim pentru Soare iar nota¸tia P o folosim pentru planete ¸si comete Metoda de calcul a coordonatelor x; y; z la un moment t, într-un sistem de coordonate iner¸tial Sxyz const¼ a în dou¼ a etape: 1) Întâi se calculeaz¼ a coordonatele orbitale polare (r; #) sau carteziene ( ; ): 2) În a doua etap¼ a se calculeaz¼ a coordonatele x; y; z Ne punem problema trecerii de la cordonatele carteziene ecliptice (x; y; z) la coordonatele carteziene ecuatoriale (x0 ; y 0 ; z 0 ): Pentru aceasta se folosesc 134 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE Figure 8 10: Orbita în sistemele de coordonate ecliptice ¸si ecuatoriale PROBLEMA LUI NEWTON 135 formulele (8 42) Este îns¼ a necesar¼ a trecerea de la elementele orbitale ecliptice ; i; ! la elementele orbitale ecuatoriale 0 ; i0 ; ! 0 (…gura 8 10) Acest lucru se face aplicând formulele lui Gauss în triunghiul sferic N N 0 : sin 0 sin 0 sin 0 sin " = sin i sin(! 0 cos " = cos(! 0 sin i 0 !) !) sin sin(! 0 !) cos cos i = sin i sin Fie X 0 ; Y 0 ; Z 0 coordonatele carteziene geocentrice ecuatoriale ale Soarelui Cunoscând coordonatele carteziene heliocentrice ecuatoriale ale astrului x0 ; y 0 ; z 0 cu formulele de transla¸tie X 0 = X 0 + x0 ; Y 0 = Y 0 + y0; Z 0 = Z 0 + z0; calcul¼ am coordonatele carteziene geocentrice ecuatoriale ale astrului În scopul compar¼ arii teoriei cu observa¸tiile, coordonatele X 0 ; Y 0 ; Z 0 se exprim¼ a cu ajutorul coordonatelor sferice geocentrice ecuatoriale R0 ; 0 ; 0 distan¸ta geocentric¼ a, ascensia dreapt¼ a geocentric¼ a ¸si declina¸tia geocentric¼ a: R0 cos 0 cos 0 = X 0 + x0 ; R0 cos 0 sin 0 = Y 0 + y0; R0 sin 0 = Z 0 + z0: Formulele precedente se folosesc pentru calculul efemeridelor planetelor, cometelor, Lunii, sateli¸tilor arti…ciali ai P¼ amântului etc În¸telegem prin efemerida unui astru un tabel cuprinzând un ¸sir de pozi¸tii ale astrului pe bolta cereasc¼ a; de exemplu un tabel cu un ¸sir de valori ale coordonatelor ecuatoriale sferice 0 ; 0 ale astrului respectiv, calculate pentru momente ce formeaz¼ a o progresie aritmetic¼ a 136 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE Capitolul 9 Problema celor n 3 corpuri Problema celor n 3 corpuri studiaz¼ a mi¸scarea a n corpuri cere¸sti P0 ; P1 ; :::; Pn de mase m0 ; m1 ; :::; mn 1 sub ac¸tiunea fortelor reciproce de atrac¸tie universal¼ a 9 1 Ecua¸ tiile mi¸ sc¼ arii absolute Fie reperul iner¸tial OXY Z; …e ri = OPi vectorul de pozi¸tie al corpului (asimilat unui punct material) Pi de coordonate (Xi ; Yi ; Zi ) ¸si …e rij = rj ri : Vom nota cu rij = rji = krj ri k = q (Xi Xj )2 + (Yi Yj )2 + (Zi Zj )2 ; (9 1) distan¸ta de la Pi la Pj Punctul Pi este atras de Pj cu for¸ta F ij = f mi mj rij : 3 rij (9 2) Ecua¸tiile de mi¸scare sub in‡uen¸ta for¸telor de atrac¸tie sunt mi n X1 rj ri d2 ri = f mi mj ; i = 0; 1; :::; n 3 2 dt rij 1: (9 3) j=0;j6=i În cazul problemei celor 3 corpuri, sistemul de ecua¸tii de mi¸scare este 137 1 138 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE d2 r0 dt2 = f m1 d2 r1 dt2 = f m0 d2 r2 dt2 = f m1 r1 r0 + f m2 3 r01 r0 r1 + f m2 3 r10 r1 r2 + f m0 3 r21 r2 r0 3 r02 r2 r1 3 r12 r0 r2 3 r20 ; (9 4) ; : În proiec¸tie pe axele sistemului de coordonate OXY Z, sistemul (9 3) are forma d2 Xi mi 2 dt = f mi mj Xj = f d2 Zi dt2 = f n X1 mi mj Yj mi mj Yi 3 rij j=0;j6=i n X1 Xi 3 rij j=0;j6=i d2 Yi mi 2 dt mi n X1 Zj j=0;j6=i (9 5) ; Zi 3 rij ; ; i = 0; 1; :::; n 1: În anul 1773, J L Lagrange (1736-1813) a introdus func¸tia de for¸t¼a a sistemului de puncte materiale: U =f n X1 nX2 j=i+1 i=0 mi mj ; rij (9 6) care în cazul n = 3 se expliciteaz¼ a astfel U =f m0 m1 m1 m2 m2 m0 +f +f : r01 r12 r20 (9 7) Sistemul (9 5) poate … pus acum sub forma d2 Xi dt2 d2 Yi mi 2 dt d2 Zi mi 2 dt mi = = = @U ; @Xi @U ; @Yi @U ; i = 0; 1; :::; n @Zi (9 8) (9 9) 1: (9 10) PROBLEMA LUI NEWTON 139 Ad¼ augând condi¸tiile ini¸tiale (0) (0) (0) Xi (t0 ) = Xi ; Yi (t0 ) = Yi ; Zi (t0 ) = Zi ; (0) (0) (9 11) (0) X i (t0 ) = X i ; Y i (t0 ) = Y i ; Z i (t0 ) = Z i ; ob¸tinem problema lui Cauchy pentru problema celor n corpuri, a c¼ arei solu¸tie exist¼ a ¸si este unic¼ a 9 2 9 2 1 Integralele prime ale sistemului de ecua¸ tii diferen¸ tiale al mi¸ sc¼ arii absolute Integralele mi¸ sc¼ arii centrului de mas¼ a Din (9 3) rezult¼ a n X1 mi n X1 mi i=0 d2 ri = 0; dt2 (9 12) de unde, prin integr¼ ari succesive ob¸tinem urm¼ atoarele integrale prime i=0 n X1 dri = A; dt mi ri = At + B: (9 13) (9 14) i=0 Proiectând rela¸tiile (9 13) ¸si (9 14) pe cele 3 axe ale sistemului de coordonate, ob¸tinem 6 integrale prime, numite integralele mi¸ sc¼arii centrului de mas¼a Din (9 14) rezult¼ a ca centrul de mas¼a al celor n corpuri, având vectorul de pozi¸tie Pn 1 mi ri = Pi=0 ; n 1 i=0 mi are o mi¸ scare rectilinie s¸i uniform¼a 9 2 2 Integralele ariilor Din (9 3) rezult¼ a de asemenea n X1 i=0 mi r i d2 ri = 0; dt2 (9 15) 140 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE ¸si integrând n X1 mi ri i=0 dri = C: dt (9 16) Proiectând (9 16) pe axele sistemului OXY Z ob¸tinem n X1 m i Yi Z i Zi Y i = CX ; m i Zi X i Xi Z i = CY ; mi Xi Y i Yi X i = CZ ; i=0 n X1 i=0 n X1 (9 17) i=0 acestea …ind cele 3 integrale prime ale conserv¼arii momentului cinetic sau integralele ariilor 9 2 3 Integrale energiei Înmul¸tind (9 8) cu X, (9 9) cu Y ; (9 10) cu Z ¸si însumând ob¸tinem n X1 mi X i X i + Y i Y i + Z i Z i i=0 = dU ; dt (9 18) de unde prin integrare ¸si cu nota¸tia n 1 E= 2 2 2 1X mi X i + Y i Z i 2 (9 19) i=0 pentru energia cinetic¼ a a sistemului, ob¸tinem integrala energiei E U = h; h …ind constanta energiei (9 20) PROBLEMA LUI NEWTON 9 3 141 Rezolvarea numeric¼ a a problemei planare a celor trei corpuri Avem de rezolvat sistemul de 6 ecua¸tii diferen¸tiale de ordinul 2: m0 d2 X0 dt2 = f m0 m1 d2 Y0 dt2 = d2 X1 dt2 = d2 Y1 dt2 = d2 X2 dt2 = d2 Y2 dt2 = m0 m1 m1 m2 m2 X1 X0 X2 X0 + f m0 m2 ; 3 3 r01 r02 Y1 Y0 Y2 Y0 f m0 m1 + f m0 m2 ; 3 3 r01 r02 X2 X1 X0 X1 f m1 m2 + f m1 m0 ; 3 3 r12 r10 Y2 Y1 Y0 Y1 f m1 m2 + f m1 m0 ; 3 3 r12 r10 X0 X2 X1 X2 f m2 m0 + f m2 m1 ; 3 3 r20 r21 Y0 Y2 Y1 Y2 f m2 m0 + f m2 m1 : 3 3 r20 r21 (9 21) Vom transforma acest sistem ïntr-un sistem de 12 ecua¸tii diferen¸tiale de 142 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE ordinul 1 : dX0 dt dY0 dt dX 0 dt dY 0 dt dX1 dt dY1 dt dX 1 dt dY 1 dt dX2 dt dY2 dt = X 0; (9 22) = Y 0; = f m1 = f m1 X1 X0 3 r01 Y1 Y0 3 r01 + f m2 + f m2 X2 X0 3 r02 Y2 Y0 3 r02 ; ; = X 1; = Y 1; = f m2 = f m2 X2 X1 3 r12 Y2 Y1 3 r12 + f m0 + f m0 X0 X1 3 r10 Y0 Y1 3 r10 ; ; = X 2; = Y 2; dX 2 dt = f m0 dY 2 dt = f m0 X0 X2 3 r20 Y0 Y2 3 r20 + f m1 + f m1 X1 X2 3 r21 Y1 Y2 3 r21 ; : Notând y = (y1 ; y2 ; y3 ; y4 ; y5 ; y6 ; y7 ; y8 ; y9 ; y10 ; y11 ; y12 )T = T X0 ; Y0 ; X 0 ; Y 0 ; X1 ; Y1 ; X 1 ; Y 1 ; X2 ; Y2 ; X 2 ; Y 2 M 0 = f m0 ; M 1 = f m1 ; M 2 = f m2 ; ; PROBLEMA LUI NEWTON 143 punem sistemul (9 22) sub forma dy1 dt dy2 dt dy3 dt dy4 dt dy5 dt dy6 dt dy7 dt dy8 dt dy9 dt dy10 dt dy11 dt dy12 dt = y3 ; (9 23) = y4 ; y5 k(y5 ; y6 ) y6 = M1 k(y5 ; y6 ) = M1 y1 (y1 ; y2 )k3 y2 (y1 ; y2 )k3 y1 k(y9 ; y10 ) y2 M2 k(y9 ; y10 ) M2 y9 ; (y1 ; y2 )k3 y10 ;; (y1 ; y2 )k3 = y7 ; = y8 ; y9 k(y9 ; y10 ) y10 = M2 k(y9 ; y10 ) = M2 y5 (y5 ; y6 )k3 y6 (y5 ; y6 )k3 y5 k(y9 ; y10 ) y6 M0 k(y9 ; y10 ) M0 y1 ; (y5 ; y6 )k3 y2 ; (y5 ; y6 )k3 = y11 ; = y12 ; y1 y9 k(y1 ; y2 ) (y9 ; y10 )k3 y2 y10 = M0 k(y1 ; y2 ) (y9 ; y10 )k3 = M0 y9 y5 ; k(y5 ; y6 ) (y9 ; y10 )k3 y10 y6 M1 : k(y5 ; y6 ) (y9 ; y10 )k3 M1 Pentru rezolvarea sistemului (9 23), vom rescrie, sub o form¼ a prescurtat¼ a, programele MATLAB prezentate de Dominik Gruntz ¸si Joerg Waldvogel în [DJ] Scriem mai întâi termenii din partea dreapt¼ a a sistemului (9 23) cu ajutorul func¸tiei simple f m: function ydot = simple f ( t, y ) %*********************************************** % Input, real T, the current time % % Input, real Y(12), the current solution % % Output, real YDOT(12), the derivatives of the current solution % 144 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE global M0 global M1 global M2 x0 = y(1:2); x1 = y(5:6); x2 = y(9:10); d0 = ( x2 - x1 ) / norm ( x2 - x1 )^3; d1 = ( x0 - x2 ) / norm ( x0 - x2 )^3; d2 = ( x1 - x0 ) / norm ( x1 - x0 )^3; ydot(1:2) = y(3:4); ydot(3:4) = M1 * d2 -M2 * d1; ydot(5:6) = y(7:8); ydot(7:8) = M2 * d0 - M0 * d2; ydot(9:10) = y(11:12); ydot(11:12) = M0 * d1 - M1 * d0; ydot = ydot(:); return end Sistemul de ecua¸tii diferen¸tiale de ordinul I (9 23) este rezolvat numeric cu metoda Runge-Kutta , folosind solverul MATLAB ode113 Pentru aceasta scriem programul three body m three body m %****************** global M0, M1, M2; % % Introduceti masele % M0 = 5 0; M1 = 3 0; M2 = 4 0; % %Fixati intervalul de timp % t initial = 0 0; t …nal = 10 0; t range = [ t initial, t …nal ]; % % Dati conditiile intiale (pozitii, viteze) PROBLEMA LUI NEWTON Figure 9 1: Solu¸tie a problemei celor trei corpuri % y initial = [ 1 0; -1 0; 0 0; 0 0; 1 0; 3 0; 0 0; 0 0; -2 0; -1 0; 0 0; 0 0 ]; % % Stabiliti eroarea maxima % options = odeset ( ’RelTol’, 1 0e-10, ’AbsTol’, 1 0E-10 ); % % Integrati sistemul de ecuatii diferentiale % [ T1, Y1 ] = ode113 (’simple f ’,t range, y initial, options ); % % A…sati rezultatele % R1 = 1 : length ( T1 ); plot ( Y1(R1,1), Y1(R1,2), ’- ’, Y1(R1,5), Y1(R1,6), ’:’, Y1(R1,9), Y1(R1,10), ’–’ ) title ( ’0 <= T <= 10’ ) 145 146 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE În …gura 9 1 prezent¼ am traiectoriile corpurilor P0 ; P1 ; P2 pentru masele ¸si conditiile men¸tionate în programul three body m De asemenea, în …gura 9 2 prezent¼ am diverse traiectorii ale corpurilor P0 ; P1 ; P2 : În …gura 9 2 a) consider¼ am masele M 0 = 10:0; M 1 = 5:0; M 2 = 0:01; ¸si condi¸tiile ini¸tiale (0) (0) (0) (0) (0) (0) (0) (0) (0) (0) (0) (0) X0 ; Y0 ; X 0 ; Y 0 ; X1 ; Y1 ; X 1 ; Y 1 ; X2 ; Y2 ; X 2 ; Y 2 = (0:0; 1:0; 1:0; 0:0; 0:0; 1:0; 2:0; 0:0; 0:0; 0:8; 0:0; 0:0) : În …gura 9 2 b) consider¼ am masele M 0 = 2:0; M 1 = 1:0; M 2 = 1:0; ¸si condi¸tiile ini¸tiale (0) (0) (0) (0) (0) (0) (0) (0) (0) (0) (0) (0) X0 ; Y0 ; X 0 ; Y 0 ; X1 ; Y1 ; X 1 ; Y 1 ; X2 ; Y2 ; X 2 ; Y 2 = (0:0; 0:0; 1:0; 0:5; 0:0; 1:0; 1:0; 0; 0:0; 1:0; 1:0; 1:0) : În …gura 9 2 c) consider¼ am masele M 0 = 2:0; M 1 = 1:0; M 2 = 1:0; ¸si condi¸tiile ini¸tiale (0) (0) (0) (0) (0) (0) (0) (0) (0) (0) (0) (0) X0 ; Y0 ; X 0 ; Y 0 ; X1 ; Y1 ; X 1 ; Y 1 ; X2 ; Y2 ; X 2 ; Y 2 = (0:0; 0:0; 0:1; 0:0; 0:0; 1:0; 1:2; 1:0; 0:0; 1:0; 1:0; 1:0) : În …gura 9 2 d) consider¼ am masele M 0 = 2:0; M 1 = 1:0; M 2 = 1:0; ¸si condi¸tiile ini¸tiale (0) (0) (0) (0) (0) (0) (0) (0) (0) (0) (0) (0) X0 ; Y0 ; X 0 ; Y 0 ; X1 ; Y1 ; X 1 ; Y 1 ; X2 ; Y2 ; X 2 ; Y 2 = (0:0; 0:0; 0:0; 0:0; 0:0; 1:0; 1:0; 1:0; 0:0; 1:0; 1:0; 1:0) : PROBLEMA LUI NEWTON 147 Figure 9 2: Problema celor trei corpuri Traiectoriile pentru diverse condi¸tii ini¸tiale ¸si diverse valori ale maselor 148 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE Capitolul 10 Elemente de trigonometrie sferic¼ a 10 1 Triunghiul sferic Triunghiul sferic ABC este …gura geometric¼ a de pe suprafa¸ta sferei , format¼ a din trei arce de cerc mari de pe sfer¼ a ce se intersecteaz¼ a ïn punctele A; B ¸si C numite vârfurile triunghiului (…gura 10 1) Arcele a; b; c;care unesc vârfurile triunghiului sferic reprezint¼ a laturile triunghiului Numim unghi al triunghiului sferic unghiul dintre tangentele la laturile sale într-unul din vârfuri, iar aceste este egal cu unghiul plan al unghiului diedru format de planele celor dou¼ a cercuri mari, a c¼ aror intersec¸tie con¸tine vârful Trei cercuri mari determin¼ a pe sfer¼ a opt triunghiuri sferice Dintre acestea com considera numai triunghiurile euleriene, adic¼ a acele triunghiuri ale c¼ aror unghiuri ¸si arce au m¼ asurile mai mici decât Între unghiurile A; B; C ¸si laturile a; b; c ale unui triunghi sferic avem urm¼ atoarele rela¸tii (pe care le d¼ am f¼ ar¼ a demonstra¸tie) c < a + b; a < b + c; b < c + a; a+b+c< ; <A+B+C <3 : Se nume¸ste exces al unui triunghi sferic m¼ arimea "= A+B+C 2 149 : 150 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE Figure 10 1: Triunghi sferic PROBLEMA LUI NEWTON 151 Figure 10 2: Triunghi sferic ¸si triunghi sferic polar 10 1 1 Triunghiul polar al unui triunghi sferic Triunghiul polar A0 B 0 C 0 al triunghiului sferic ABC este un triunghi sferic ale c¼ arui vârfuri sunt polii sferici ai cercurilor mari din care fac parte arcele BC; CA ¸si AB (…gura 10 2) Avem deci OA0 ? (OBC) ; OB 0 ? (OAC) ; OC 0 ? (OAB) (10 1) OA ? OB 0 C 0 ; OB ? OA0 C 0 ; OC ? OA0 B 0 : (10 2) ¸si de aici A¸sadar triunghiul polar al triunghiului polar este triunghiul sferic ini¸tial Între elementele triunghiului sferic ABC ¸si elementele triunghiului s¼ au polar A0 B 0 C 0 avem rela¸tiile A + a0 = B + b0 = C + c0 = 1800 ; (10 3) A0 + a = B 0 + b = C 0 + c = 1800 : (10 4) 152 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE Figure 10 3: Pentru deducerea formulelor lui Gauss 10 2 Formulele lui Gauss Fie (vezi …gura 10 3) triunghiul sferic ABC pe sfera de centru O Vom considera dou¼ a sisteme de coordonate carteziene cu originea în O Primul sistem Oxyz are axa Oz în direc¸tia OA Consider¼ am Oy în plaul OAB astfel încât Oy ? Oz iar Ox ? yOz, de aceea¸si parte cu C Al doilea sistem Ox0 y 0 z 0 are axa Oz 0 în direc¸tia OB Consider¼ am Oy 0 în plaul OAB astfel 0 0 0 0 încât Oy ? Oz iar Ox Ox Fie D ¸si D proiec¸tiile lui C pe planele xOy 0 0 respectiv x Oy :Avem rela¸tiile \0 sin COB; [ sin COA; [ x0 = kOCk cos xOD \ (10 5) x = kOCk cos xOD \0 sin COB; [ sin COA; [ y 0 = kOCk sin xOD \ y = kOCk sin xOD [ z 0 = kOCk cos COB: \ z = kOCk cos COA; Dar [ =A xOD 2 \0 = ; xOD 2 [ = b; COB \ = a: B; COA (10 6) Cum Oy; Oy 0 ; OB ¸si OA sunt coplanare ¸si Oy ? OA; Oy 0 ? OB, deducem PROBLEMA LUI NEWTON 153 c¼ a [0 = BOA \ = c: yOy A¸sadar, pentru a suprapune Oy 0 peste Oy ¸si Oz 0 peste Oz trebuie s¼ a facem o rota¸tie de unghi ( c) în jurul axei Ox: Avem deci x0 = x; (10 7) y 0 = y cos c + z sin c; z 0 = y sin c + z cos c: Înlocuind x; y; z; x0 ; y 0 ; z 0 în (10 7) cu expresiilor din (10 5) ¸si ¸tinând cont de (10 6) ob¸tinem formulele lui Gauss: sin a sin B = sin b sin A; sin a cos B = cos b sin c (10 8) sin b cos c cos A; cos a = cos b cos c + sin b sin c cos A; 10 2 1 Formulele sinusurilor Prin permut¼ ari circulare, prima din formulele (10 8) conduce la sin a sin B = sin b sin A; (10 9) sin b sin C = sin c sin B; sin c sin A = sin a sin C: 10 2 2 Formulele care cuprind cinci elemente A doua dintre rela¸tiile (10 8), la care se adaug¼ a înc¼ a cinci rela¸tii asem¼ an¼ atoare, ob¸tinute prin permut¼ ari ale elementelor triunghiului sferic, ¸tinând cont c¼ a …ecare latur¼ a a triunghiului poate … asociat¼ a, pe rând, cu câte unul din cele dou¼ a unghiuri adiacente, constituie formulele care cuprind cinci elemente sin a cos B = cos b sin c sin b cos c cos A; sin a cos B = cos c sin b sin c cos b cos A; sin b cos C = cos c sin a sin c cos a cos B; sin b cos C = cos a sin c sin a cos c cos B; sin c cos A = cos a sin b sin a cos b cos C; sin c cos A = cos b sin a sin b cos a cos C; (10 10) 154 10 2 3 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE Formulele cosinusurilor Prin permut¼ ari circulare, prima din formulele (10 8) conduce la cos a = cos b cos c + sin b sin c cos A; (10 11) cos b = cos c cos a + sin c sin a cos B; cos c = cos a cos b + sin a sin b cos C: Rescriind formulele (10 11) pentru elementele triunghiului polar A0 B 0 C 0 avem cos a0 = cos b0 cos c0 + sin b0 sin c0 cos A0 ; 0 0 cos b = cos c cos a + sin c sin a cos B ; 0 = cos a0 cos b0 + sin a0 sin b0 cos C 0 : cos c 0 0 0 (10 12) 0 Aplicând formulele (10 3) ¸si (10 4) în (10 12) ob¸tinem formulele cosinusurilor pentru unghiurile triunghiului sferic cos A = cos B cos C + sin B sin C cos a; cos B = cos C cos A + sin C sin A cos b; 0 cos A cos B + sin A sin B cos c: cos C 10 2 4 = (10 13) Expresiile unghiurilor în func¸ tie de laturi Din prima formul¼ a a cosinusurilor ob¸tinem cos A = a+b c a b+c sin A 2 2 = 2 sin2 ; sin b sin c 2 2p = a + b + c; 2 sin 1 cos A = A sin = 2 2 sin 1 + cos A = cos a cos b cos c ; sin b sin c r sin (p b) sin (p sin b sin c a+b+c b+c sin 2 2 sin b sin c c) ; (10 14) a = 2 cos2 A ; 2 PROBLEMA LUI NEWTON r A sin p sin (p a) cos = ; 2 sin b sin c s sin (p b) sin (p c) A : tan = 2 sin p sin(p a) Asem¼ an¼ ator, folosind celelalte formule ale cosinusurilor rezult¼ a r B sin (p c) sin (p a) sin = ; 2 sin c sin a r sin p sin (p b) B ; cos = 2 sin c sin a s B sin (p c) sin (p a) tan = ; 2 sin p sin(p b) r sin (p a) sin (p b) C ; sin = 2 sin a sin b r C sin p sin (p c) cos = ; 2 sin a sin b s sin (p a) sin (p b) C : tan = 2 sin p sin(p c) 10 2 5 155 (10 15) (10 16) (10 17) (10 18) (10 19) (10 20) (10 21) (10 22) Expresiile laturilor în func¸ tie de unghiuri Din prima formul¼ a a cosinusurilor pentru unghiurile triunghiului sferic se ob¸tine cos A + cos B cos C ; sin B sin C A+B+C B+C A 2 cos cos a 2 2 1 cos a = = 2 sin2 ; sin B sin C 2 2" = A + B + C ; r sin " sin(A ") a sin = ; 2 sin B sin C B C +A B C A 2 cos cos a 2 2 1 + cos a = = 2 cos2 ; sin B sin C 2 cos a = (10 23) 156 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE a cos = 2 tan a = 2 r s sin(B ") sin(C sin B sin C ") ; (10 24) sin " sin(A ") : sin(B ") sin(C ") (10 25) Cu celelalte dou¼ a formule se ob¸tine r b sin " sin(B ") ; sin = 2 sin C sin A r b sin(C ") sin(A ") cos = ; 2 sin C sin A s b sin " sin(B ") tan = ; 2 sin(C ") sin(A ") r c sin " sin(C ") sin = ; 2 sin A sin B r sin(A ") sin(B ") c ; cos = 2 sin A sin B s sin " sin(C ") c tan = : 2 sin(A ") sin(B ") 10 2 6 Formulele lui Delambre Din formulele (10 14)-(10 30), ob¸tinem formulele lui Delambre cos A B 2 C sin 2 sin A 2 C cos 2 a+b 2 ; c sin 2 sin = B sin = a 2 c sin 2 b ; A+B a+b cos 2 2 ; = c C cos sin 2 2 cos (10 26) (10 27) (10 28) (10 29) (10 30) (10 31) PROBLEMA LUI NEWTON 157 A+B a b cos 2 2 : = c C cos cos 2 2 sin 10 2 7 Formula lui L’Huillier Vom exprima excesul triunghiului sferic în func¸tie de laturi Pornim de la formula x y sin x sin y tan = : 2 cos x + cos y Cum " A+B+C tan = tan 2 4 avem cu formula de mai sus A+B sin sin " 2 tan = A+B 2 cos + cos 2 A+B 2 = tan C C 2 2 ; A+B C cos 2 2 2: = C A+B C cos + sin 2 2 2 sin A+B A+B Înlocuind sin ¸si cos cu formulele corespunz¼ atoare ale lui De2 2 lambre rezult¼ a a b c cos cos " 2 2 cot C ; tan = a + b c 2 2 cos + cos 2 2 sau a b+c c a+b sin sin " C 4 4 tan = cot ; a+b+c a+b c 2 2 cos cos 4 4 deci cu formula (10 31) p a p bs sin sin " sin p sin (p c) 2 2 tan = : p p c 2 sin (p a) sin (p b) cos cos 2 2 Ridicând la p¼ atrat ob¸tinem formula lui Hillier tan2 " p p a p b p c = tan tan tan tan : 2 2 2 2 2 158 DINAMICA PUNCTULUI S¸I A SISTEMELOR DE PUNCTE 10 2 8 Aria triunghiului sferic Fie A0 ; B 0 ; C 0 punctele diametral opuse punctelor A; B; C de pe sfer¼ a (…gura ??) Avem rela¸tiile: b SABC + SBCA0 = 2R2 A; b (aria fusului sferic de unghi A), b SABC + SCAB 0 = 2R2 B; b (aria fusului sferic de unghi B), b SABC + SABC 0 = 2R2 C; b (aria fusului sferic de unghi C) Adunând formulele de mai sus, ob¸tinem b+B b+C b : 3SABC + SBCA0 + SCAB 0 + SABC 0 = 2R2 A (10 32) Datorit¼ a simetriei fa¸ta¼ de centru, avem SABC + SBCA0 + SCAB 0 + SABC 0 = SA0 B 0 C 0 + SB 0 C 0 A + SC 0 A0 B + SA0 B 0 C : Triunghiurile din rela¸tia precedent¼ a acoper¼ a întreaga suprafa¸ta¼ a sferei SABC +SBCA0 +SCAB 0 +SABC 0 +SA0 B 0 C 0 +SB 0 C 0 A +SC 0 A0 B +SA0 B 0 C = 4 R2 ; deci SABC + SBCA0 + SCAB 0 + SABC 0 = 2 R2 : Sc¼ azând (10 33) din (10 32) rezult¼ a deci b+B b+C b 2SABC = 2R2 A b+B b+C b SABC = R2 A unde " este excesul triunghiului sferic ; = 2R2 "; (10 33) Bibliogra…e [DJ] Dominik Gruntz, Joerg Waldvogel, "Orbits in the Planar Three-Body Problem",in Walter Gander, Jiri Hrebicek, Solving Problems in Scienti…c Computing using Maple and Matlab, Springer, 1997, ISBN: 3-54061793-0 159 